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INTRODUCTION 


The  analytical  determination  of  the  pressure  acting  on  the  blunt 
base  of  axisymmetric  supersonic  vehicles  has  been  under  consideration 
for  more  than  two  decades.  Although  a  successful  analytical  approach 
for  the  plane  two-dimensional  turbulent  supersonic  base  pressure  problem 
with  uniform  external  flow  approaching  the  base  evolved  in  the  early 
1950* s,  an  equally  successful  and  straightforward  analytical  method  for 
the  axisymmetric  case  was  not  developed  until  the  late  1960's.  In  a 
review  of  base  drag  research  by  Sedney1  in  1966,  it  was  concluded  that 
no  satisfactory  theory  existed.  The  principal  objectives  of  the  work 
described  herein  were  to  modify  and  improve  the  method  of  calculating 
the  turbulent  base  pressure  in  axisymmetric  supersonic  flow,  develop  a 
computer  program  using  this  method,  and  compare  computed  results  with 
those  from  experiments  for  cases  of  interest  to  the  U.S.  Army  Ballistic 
Research  Laboratory,  ARRADCOM. 


AXISYMMETRIC  BASE  PRESSURE  PROBLEM 

Analytical  studies  were  conducted  to  evaluate  the  utility,  for  the 
present  task,  of  the  mathematical  base  pressure  model  for  axisymmetric 
supersonic  flow  as  developed  by  Zumwalt*  and  modified  by  Mueller, 
et  al3’4.  It  was  found  that  this  turbulent  flow  model  offered  the  best 
possibility  of  predicting  the  base  pressure  for  high  Reynolds  number 
flows  and  for  including  the  effects  of  base  bleed  and  Reynolds  number 
later  with  a  minimum  of  computer  time.  The  configurations  of  interest 
in  this  investigation  are  shown  in  Figure  1.  For  high  Reynolds  numbers, 
the  turbulent  boundary  layers  on  the  short  bodies  shown  in  Figure  1  are 
usually  very  thin,  while  the  inviscid  flow  upstream  of  the  base  in  la 
and  Id  or  upstream  of  the  afterbody  in  lb  and  lc  may  be  uniform  or 
non-uniform. 


1.  Sedney,  R. ,  "Review  of  Base  Drag",  BRL  Report  No  1337,  U.S.  Amy 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland, 
October  1966.  AD  808767 

2.  Zumwalt,  G.  W. ,  "Analytical  and  Experimental  Study  of  the 
Axially-syrmetric  Supersonic  Base  Pressure  Problem,  Ph.D. 
Dissertation" ,  University  of  Illinois,  1959. 

3.  Mueller,  T.  J.,  "Determination  of  the  Turbulent  Base  Pressure  in 
Supersonic  Axisymmetric  Flow".  Journal  of  Spacecraft  and  Rockets, 
Vol.  5,  No.  1,  pp.  101-107,  January  1968. 

4.  Mueller,  T.  J.,  and  Hall,  C.  R.,  Jr.,  "Analytical  Prediction  of  the 
Turbulent  Base  Pressure  in  Supersonic  Axisymmetric  Flow  Including 
the  Effect  of  Initial  Flow  Direction".  AFFDL-TR-68-132, 

September  1968. 
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The  computer  program,  called  BOATTAIL/ FLARE -BASE  5,  uses  the 
axisymmetric  method  of  characteristics  to  calculate  the  flow  over 
conical  boattailed  and  flared  bodies  of  revolution.  This  calculation 
yields  the  necessary  initial  flow  field  for  the  base  pressure  determi¬ 
nation.  This  program  can  be  used  to  calculate  the  turbulent  base 
pressure  for  configurations  similar  to  all  those  in  Figure  1.  Both 
non-uniform  and  uniform  supersonic  flows  approaching  the  base  or 
afterbody  can  be  handled. 

The  theoretical  flow  model  used  in  conjunction  with  the  axisymmetric 
method  of  characteristics  for  the  flow  over  boattails  and  flares  is 
shown  in  Figure  2.  In  this  model,  the  flow  is  divided  into  three 
components:  (a)  an  inviscid  free  stream,  (b)  a  dissipative  mixing  layer 
and  (c)  a  base  region.  In  addition  to  using  the  restricted  mixing 
theory  of  Korst5,  the  following  conditions  are  imposed  on  these  three 
flow  components: 

(i)  The  boundary  layer  approaching  the  separation  corner  is 
neglected  although  it  is  assumed  to  be  fully  turbulent. 

(ii)  An  issentropic  expansion  takes  place  in  the  free  stream 
from  1  to  2,  i.e.,  a  continuation  of  the  axisymmetric 
method  of  characteristics  for  the  flow  over  the  afterbody. 

(iii)  The  inviscid  flow  past  a  conetail  using  axisymmetric 
characteristics  is  utilized  to  define  the  pressure  field 
impressed  on  the  mixing  region  from  2  to  3.  The  conetail 
surface  also  serves  as  the  "corresponding  inviscid  jet 
boundary" . 

(iv)  The  pressure  normal  to  the  "corresponding  inviscid  jet 
boundary"  is  assumed  to  be  constant  within  and  near  the 
mixing  region  at  each  cross-section. 

(v)  Velocity  profile  similarity  is  assumed  in  the  mixing 
region.  The  error  function  velocity  distribution  is 
located  within  the  intrinsic  system  of  coordinates  x,y. 

The  coordinate  shift  between  the  reference  coordinates 
X, Y  and  the  intrinsic  coordinates  x,y  is  represented  by 
X  -  x  and  Y  =  y-ym(x)  where  ym(0)  =  0.  This  coordinate 

shift  is  a  result  of  using  the  restricted  mixing  theory 
of  Korst. 


5.  Korst ,  H.H.,  "A  Theory  for  Base  Pressures  in  Transonic  and 

Supersonic  Flow".  Journal  of  Applied  Mechanics,  Vol.  23,  No.  4, 
pp.  593-600,  December  1956. 


(vi)  The  axisymmetric  geometry  of  the  mixing  region  is  taken 
into  account  in  the  integral  presentation  for  momentum  and 

mass  flux  between  sections  2  and  3.  j 

(vii)  Recompression  is  assumed  to  result  from  an  oblique  shock  . 

turn  from  3  to  4  at  the  empirically  determined  trailing 

wake  radius.  Values  for  the  flow  over  cylindrical  bodies 
are  given  by  Chapman6,  and  are  reproduced  in  Figure  3. 


A  streamline,  j,  can  be  identified  which  divides  the  amount  of  mass 
passing  over  the  corner  at  1  from  that  mass  flow  entrained  by  the 
viscous  action  of  the  free  jet  mixing  region.  A  second  streamline,  d, 
can  be  identified  which  has  just  sufficient  kinetic  energy  at  3  to 
negotiate  the  pressure  rise  to  4.  Streamlines  below  the  d-streamline 
have  lower  kinetic  energies  and  are  unable  to  pass  through  recompression 
and  are  turned  back  to  recirculate  in  the  base  region.  If  there  is  no 
secondary  flow,  i.e.,  no  mass  bleed  into  or  out  of  the  base  region,  the 
conservation  of  mass  requires  that  the  j-  and  d-streamlines  be  identical. 
Korst  distinguished  between  these  two  streamlines  and  pictured  the  space 
between  them  as  a  sort  of  corridor  through  which  mass  would  flow  into 
or  out  of  the  base  region. 


The  control  volume  utilized  in  the  formulation  of  the  conservation 
of  mass  and  momentum  between  cross-sections  2  and  3  is  bounded  by 
streamlines  R  and  -R  as  shown  in  Figure  2.  The  streamlines  R  and  -R 
were  defined  by  Zumwalt2  so  that  the  cross-sectional  area  normal  to  the 
direction  of  flow  would  remain  nearly  constant  and  the  P  dA  pressure 
force  could  be  neglected  in  the  momentum  equation.  This  condition  is 
retained  in  the  present  development  although  its  validity  is  question¬ 
able.  For  the  simplified  axisymmetric  flow  model  described  above, 
Zumwalt  formulated  the  momentum  equation  in  the  axial  direction  using 
geometrical  relations  and  the  relation  between  the  inviscid  and  viscid 
coordinate  systems.  This  equation  was  solved  simultaneously  with  the 
combined  viscid  and  inviscid  continuity  equations  written  for  the 
control  volume  between  cross-sections  2  and  3.  For  the  error  function 
velocity  profile,  =  1/2(1  +  erf  n)  where  4  =  u/ug  and  n  =  oy/x,  it 

was  found  that  nR  =  3  was  large  enough  for  to  approach  its  asymptotic 


value.  The  result  of  this  analysis  is  a  nonlinear  equation  which  allows 
one  to  locate  the  j -streamline  at  cross-section  3,  namely: 


6.  Chapman ,  D.R.,  "An  Analysis  of  Base  Pressure  at  Supersonic 

Velocities  and  Comparison  with  Experiment".  NACA  Report  1051,  1951. 
(Supercedes  NACA  TN-2137,  1950.) 
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It  is  important  to  notice  that,  since  these  integrals  are  expressed 
only  in  terms  of  Crocco  number,  they  are  independent  of  y  although 
Equation  1  is  a  function  of  y. 

It  is  evident  that,  to  determine  4> ^ 3  from  Equations  1  and  2  for  a 
given  initial  condition,  the  location  of  the  reconroression  revion, 
r3/r^,  the  corresponding  inviscid  condition,  M3a  and  the  jet  spread 
parameter,  o3a,  must  be  known.  The  location  of  the  recompression 
region,  r3/rb>  is  determined  from  the  experimental  data  for  cylindrical 

bodies  and  is  replotted  in  Figure  3.  Unfortunately,  the  effect  of 
initial  flow  direction  on  the  wake  radius  ratio  does  not  appear  to 
have  been  documented  experimentally  or,  if  it  has,  it  has  not  been 
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published  in  the  open  literature, 
jet  boundary  from  2  to  3,  i.e.,  M 


The  Mach  number  along  the  inviscid 


3a’ 


is  determined  from  the  axisymmetric 


characteristics  solution.  The  jet  spread  parameter,  03a>  is  determined 
from  the  equation  obtained  by  Channapragada6: 


o3a  =  { R[ 1  ♦  8(1  -  C3  ) ] ) 


'  inc 


(71 


where  a.  =12  for  the  error  function  velocity  profile  used, 
inc  r 

B  =  TQa/Tjj  =  1  for  isoenergetic  mixing  considered  here,  and  R  is  the 

empirical  compressible  divergence  factor  defined  by  Channapragada7  and 
presented  as  a  function  of  Crocco  number. 

2 

The  geometric  parameter,  [or/x  cos8 ] 3  ,  may  be  calculated,  since 

a 

03_4  =  0i-2  from  the  conical  wake  assumption  and  since 


or 

2 

0  tare 

x  cose. 

3a 

L[l/(f/rb)]-l 

(8) 


At  this  point  4j3  =  <}>d3 


(i.e.,  no  base  bleed)  may  be  obtained  from 
Equations  1  and  2;  therefore,  Cd3  =  ^dsSa  =  j 3C3a  is  calculated  for 
isoenergetic  flow.  However,  the  value  of  Cd3  can  also  be  obtained  from 

the  assumed  isentropic  recompression  mechanism  along  the  d-streamline 
(i-e. ,  P 


o3d 


P4) ,  from 


Cd3  =  U  *  Cp4/pj)  'l(Y  ‘  1)/YV/2 


(9) 


where  P4/P3  is  the  pressure  rise  across  the  two-dimensional  oblique 

shock  which  results  from  the  turning  of  the  flow  with  Mach  number  MJa 

through  the  angle  03_i When  the  two  values  of  Cd3  are  equal,  then 

the  assumed  base  pressure  ratio  is  the  correct  one  for  the  given 
conditions  at  1. 

The  flow  model  used  for  constant  pressure  jet  mixing  is  identical 
to  the  flow  model  used  in  the  previous  section,  except  that  assumption 
(iii)  of  the  previous  section  is  discarded  and  replaced  by  the 


7,  Channapragada,  R.S.,  "Compressible  Jet  Spread  Farameter  for 

Mixing  Zone  Analyses".  AIAA  Journal,  Vol.  1,  No.  9,  pp.  2188-2189, 
September  1963. 
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assumption  that  the  velocity  of  the  jet  adjacent  to  (but  outside  of) 
the  mixing  region  is  constant  in  magnitude  although  not  in  direction 
For  this  situation,  61-2  *63-4  and  the  shape  of  the  "corresponding 
inviscid  jet  boundary"  is  determined  from  the  axisymmetric  method  of 


is  obtained  from  Equations  1  and  2  by  setting  Pfe=  P2=Pj  an<*  C3a=C2a' 

FORTRAN  IV  COMPUTER  PROGRAM 

A  FORTRAN  IV  computer  program  for  the  UNIVAC  1108  digital  computer 
was  developed  to  calculate  the  axisymmetric  flow  field  over  conical 
boattailed  and  flared  bodies  of  revolution  using  the  axisymmetric 
method  of  characteristics.  The  results  of  this  calculation  provide 
the  necessary  initial  flow  field  for  the  base  pressure  determination 
which  forms  the  second  part  of  the  total  computer  program.  Complete 
program  symbols,  usage  and  listing  are  presented  in  Appendices  A,  B 
and  C. 


COMPUTER  PROGRAM  DESCRIPTION 

Because  of  the  complex  and  unique  nature  of  this  program,  a 
serious  effort  directed  toward  understanding  its  operation  and 
limitations  is  recommended  before  it  is  run  on  the  computer. 

Method  of  Characteristics  for  Boattails  and  Flares 


A  schematic  of  an  axisymmetric  characteristic  grid  is  shown  in 
Figure  4.  Values  of  the  variables  X,R,v,  and  0  are  determined  at  the 
1=1  position,  and  the  calculation  of  the  remainder  of  the  points  is 
obtained  from  these  initial  conditions  and  the  boundary  conditions 
along  the  walls.  The  calculation  procedure  sweeps  diagonally  as  shown 
by  the  arrow  in  Figure  4.  The  calculations  begin  at  a  point  off  the 
surface  of  the  body  revolution,  and  proceed  outward  until  the  limit  of 
the  matrix  is  reached  at  the  1=1  location.  The  limit  referred  to  here 
is  the  limit  on  J  since  I  is  a  constant  (  =  1)  at  this  location.  For 
the  interior  flow  field  points  (i.e.,  those  points  not  on  the  surface 
of  the  body),  the  subroutine  KBPR.CALC  is  called  which  determines  the 
values  of  the  variables  at  a  new  grid  point  further  downstream.  From 
Figure  5a  it  is  noted  that  for  I  an  odd  number,  the  position  of  the 
point  (I+1,J)  may  be  determined  from  the  left  running  wave  of  the 
(I , J - 1 )  point  and  the  right  running  wave  of  the  (I,J)  grid  point.  Once 
the  X  and  R  coordinates  of  this  new  (I+1,J)  grid  point  have  been 
established,  the  additional  variables  of  v  and  0  may  be  determined. 

For  I  being  an  even  number,  the  (I,J)  and  (I,J-1)  grid  points  are  used 
to  determine  the  (I+l.J-l)  grid  point  (see  Figure  5b). 
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The  key  to  the  mesh  system  of  the  characteristic  solution  is  that 
there  are  no  J=1  points  when  I  is  an  even  number  (see  Figure  6).  It 
should  be  noted  that  the  first  number  in  the  parenthesis  will  be  the 
value  of  the  subscript  I,  while  the  second  number  refers  to  the 
subscript  J.  In  doubly-subscripted  variables,  I  is  along  the  axis  of 
the  body,  while  J  is  in  the  radial  direction.  For  singly-subscripted 
variables,  the  subscript  I  is  not  confined  to  any  particular  location, 
but  instead  refers  to  the  order  of  points. 

On  the  body  surface  and  on  the  conetail  surface  in  the  base 
pressure  solution,  the  location  of  a  characteristic  grid  point  is 
obtained  by  use  of  the  subroutine  KBPR. SURF.  Since  the  J=1  point 
always  falls  on  the  surface  itself,  only  the  (1-1,2)  point  is  required 
to  locate  the  new  (1,1)  point.  This  is  because  the  value  of  8  at  the 
surface  is  known  (i.e.,  8  at  the  surface  is  parallel  to  the  surface). 

The  characteristic  solution  continues  until  a  point  on  the  surface 
exceeds  the  end  point  of  the  body  (or  the  end  point  of  the  portion  of 
the  body  being  considered),  XFIN.  From  this  point,  the  subroutine 
KBPR. SETUP  is  called  in  an  effort  to  locate  a  grid  point  on  the  body  as 
close  to  this  end  point  as  possible.  Actually,  KBPR. SETUP  is  called 
twice  to  further  refine  the  initial  correction.  Once  the  characteristic 
line  to  this  point  has  been  established,  the  characteristics  solution 
proceeds  to  line  (I)  or  line  (2)  as  shown  in  Figure  4.  Location  (l) 

(or  Q)  )  is  read  into  the  program  depending  on  the  type  of  output 
desired.  Location  (I)  refers  to  a  constant  I  output  (not  necessarily 
a  straight  line),  while  location  (?)  is  the  output  along  a  characteristic 
line.  Again,  it  should  be  emphasized,  that  even  though  the  last  point 
on  the  surface,  XFIN,  has  been  determined,  the  characteristics 
solution  still  proceeds  until  the  values  of  J  have  been  exhausted  at 
the  1=1  position.  The  characteristics  solution  then  proceeds  only  to 
this  line  (!)  (or  Q)  ),  thus  eliminating  needless  calculations. 

Once  the  characteristics  solution  has  been  completed  for  the 
section  of  the  body  of  revolution  being  considered,  the  values  of  the 
variables  X,R,v  and  8  along  (l)  (or  (2)  )  are  tabulated  separately.  To 
avoid  confusion,  these  variables  are  now  referred  to  as  X2,  R2,  v2  and 
92  (see  Figure  7).  When  XFIN  is  greater  than  the  starting  value  of 
X ( 1 , 1 ) ,  there  are  always  fewer  values  of  the  variables  at  any  location 
in  the  characteristics  matrix  than  there  are  at  the  1=1  location.  This 
is  because  of  the  inherent  nature  of  the  characteristics  solution  at  a 
free  boundary  (see  Figure  4).  The  tabulated  values  of  X2,  R2,  v2  and 
92  are  now  stored  and  may  be  used  to  further  proceed  with  another 
characteristics  solution  for  another  downstream  section  of  the  body  of 
revolution.  It  should  be  noted  that  these  variables  need  only  one 
subscript  since  the  line  (I)  (or  (2)  )  describes  the  location  in  space. 
The  subscript  merely  refers  to  the  order  of  the  points,  i.e.,  1=1 
refers  to  the  point  on  the  body,  1=2  is  the  next  point  out,  etc. 
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An  independent  variable  is  now  introduced,  and  the  values  of  X2, 
R2,  v2  and  02  are  assumed  to  be  functions  only  of  this  new  variable,  D. 
More  precisely,  D  measures  the  distance  from  the  body  surface  along  the 
line  (T)  (or  (T)  ).  This  independent  variable,  D,  is  defined  as: 

D(I)  =  D(I-l)  ♦  /[X2(I)  -  X2 ( I  —  1 ) ] ^  +  [R2(I)  -  R2(I-1)]^ 


This  equation  shows  that  D(I)  is  the  linear  distance  of  some  grid  point 
from  the  preceding  grid  point  plus  the  sum  of  the  linear  distances  from 
previous  grid  points.  Thus  the  following  relations  of  X2,  R2,  v2  and 
02  vs.  D  are  shown  qualitatively  in  Figure  8  for  the  case  of  the  output 
being  along  a  characteristic  (line  @  ).  Therefore,  for  any  distance 
along  (2)  (or  (l)  )  there  are  values  of  X2,  R2,  v2  and  02  available. 

The  characteristics  solution  now  proceeds  to  the  next  section  of 
the  body  of  revolution.  Before  the  solution  continues,  however,  a  full 
starting  matrix  of  X,R,v  and  0  must  be  determined  at  the  new  1=1 
position.  This  is  obtained  by  defining  a  new  independent  variable, 

D1(I).  Both  D1 (I)  and  D(I)  are  distances  measured  along  (I)  (or  (5)  ). 
However,  D1(I)  defines  the  grid  spacing  required  to  re-start  the 
characteristics  solution  for  the  new  region.  These  new  grid  spacings 
are  functions  of  the  boattail  (or  flare)  angle,  and  also  the  surface 
Mach  number.  Once  these  values  of  D1  have  been  determined,  values  of 
X,R,v  and  0  are  obtained  along  (l)  (or  (g)  )  by  interpolating  the 
values  of  X2,  R2,  v 2  and  0 2  of  Figure  8.  The  KBPR.TAB  sub-program 
which  allows  any  order  interpolation  up  to  order  10  is  used  to  carry 
out  this  operation.  The  characteristics  solution  is  now  set  to  proceed 
through  this  new  region  of  the  body.  The  calculation  procedure  is 
identical  to  that  described  above,  and  only  minor  changes  in  computer 
logic  are  made  with  each  section  (except  the  base  pressure  solution). 

This  above  discussion  of  the  characteristics  solution  refers 
primarily  to  the  subroutine  KBPR. FLARE  which  directs  the  actual 
calculation  of  the  characteristics  solution  along  the  body  of  revolution. 
Flow  Chart  1  in  Appendix  D  indicates  the  operation  of  KBPR. FLARE. 

Base  Pressure  Solution 

The  calculation  procedure  in  the  base  pressure  section  of  the 
computer  program  may  be  divided  into  seven  major  steps.  The  following 
description  of  the  calculation  procedure  is  taken  essentially  from 
reference  2.  The  flow  model  is  shown  in  Figure  2. 

Step  1:  A  base  pressure  ratio  P^/P ^  (  =  PBP1)  is 
assumed  and  the  corresponding  value  of  the 
inviscid  jet  boundary  angle  (THET12) ,  and 
M2A  are  obtained  from  a  Prandtl -Meyer 
expansion  from  region  I  to  2.  The  Crocco 
number,  C2A,  is  obtained  from  the  relation: 
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(10) 


C*A 


m2a 


Step  2:  The  value  of  the  wake  radius  ratio,  R^/R^  (=  RSRB) ,  is 

obtained  from  the  function  subprogram  KBPR.RSBF.  The 

Mach  number  ratio  M3  /M2,  (=  M3AM2A)  is  obtained  from  an 
Si  a 

axisymmetric  method  of  characteristics  solution  of  the 
near  wake  flow  field.  Then  the  Crocco  numbers,  C3A  and 
C3A/C2A,  are  determined  as  functions  of  M3A  and  M2A;  also 
the  angle  of  the  inviscid  jet  boundary,  THET3A,  is 
obtained. 

Step  3:  The  jet  spread  parameter  at  the  recompression,  S3A,  is 
calculated  (in  the  subprogram  KBPR. SIGMA)  as  a  function 
of  M3A.  The  geometric  parameter,  GP,  is  then  calculated 
using : 


S3A  x  tan(THET12)1  2 

,  (1/RSRB)  -  1  J  UiJ 

Step  4:  The  velocity  ratio  PHIJ3  is  now  solved  from  a  complicated 
functional  relationship: 

PHIJ3  «  <».  ■  fen  (GP,  M3A,  G,  C3A/C2A) 

3 J 

Specifically,  <)> .  is  found  from: 

J  3 

*j3  =  I  (l  +  erf  ^j,)  (12) 

where  is  the  value  of  the  upper  integral  limit 
which  satisfied: 


KBPR. FSTEP4 


( (B-3) 2  +  2 (1-C2) [BI j 


-  Jl 


’J 


])  -  GP  =  0 


where  B  is  given  by  Equation  2,  and  the 
integrals  I x ,  I2,  Ji  and  J2  are  given  by 


Equations  3,  4,  S  and  6,  respectively. 
The  error  function  may  be  expressed: 


erf (n) 


and  erf  (n)  is  evaluated  in  the  function 
subprogram  KBPR. ERF. 


(13) 


(14) 
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The  pressure  ratio,  P2/P3,  is  determined 
from  isentropic  relations,  and  from  the  Mach 
number  ratio,  M2A/M3A,  obtained  in  Step  2. 

Step  5:  For  steady  flow  with  no  mass  bleed,  Cj3  (=C5) 
may  be  determined  from: 


Cd3  "  CS  =  $j3C3A 


(15) 


Step  6:  The  pressure  ratio,  P4/P3,  is  then  calculated 

from  the  oblique  shock  relations  with  the  initial 
Mach  number  being  M3A  and  the  turning  angle  being 
THET3A.  With  isentropic  recompression  assumed 
along  the  dividing  streamline,  then: 

1 


“d3 


e  C6  = 


1  - 


1  1 2 


(P4/P3) 


Hi 

Y 


(16) 


Step  7:  The  subroutine  KBPR. BASES  takes  the  difference 
between  C^3  s  C5  (obtained  in  Step  5)  and 

Cj3  e  C6  (obtained  in  Step  6).  When  CS  =  C6, 

a  solution  is  obtained.  New  estimates  of  the 
base  pressure  are  made  by  using  the  subroutine 
KBPR. LINEAR  until  CS  =  C6,  indicating  a 
solution  has  been  reached.  The  basic  convergence 
scheme  for  this  base  pressure  analysis  is  shown 
schematically  in  Flow  Chart  3  in  Appendix  D. 

A  general  Flow  Chart  of  the  base  pressure  section 
of  the  report  is  shown  in  Flow  Chart  4  in  Appendix  D. 

There  are  two  possible  solutions  which  may 
be  used  in  the  base  pressure  analysis.  Both 
involve  axisymmetric  method  of  characteristics 
solutions.  These  options  are: 

a)  Conetail  solution  (IBOUND-1)  in  which 
THET12  =  THET3A.  However,  M3 A  *  M2A. 

b)  Constant  pressure  boundary  solution 
(I BOUND® 2)  in  which  M3A  =  M2A.  In 
this  case,  however,  THET12  *  THET3A. 

In  some  cases,  the  stability  of  the 
program  depends  on  the  initial  estimate  of  the 
base  pressure,  PEST,  which  is  read  into  the 
program.  In  general,  a  flared  body  requires 
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a  more  accurate  guess  than  a  boattailed  after¬ 
body.  If  the  range  of  a  solution  is  not  known, 
it  is  better  to  estimate  a  low  base  pressure 
ratio  for  flared  bodies.  For  boattailed 
afterbodies,  much  more  latitude  is  permissible 
in  estimating  the  base  pressure  ratio.  However, 
for  large  boattail  angles  (>  12°),  a  high 
estimate  of  PEST  should  be  used. 

High  Mach  numbers  also  present  some 
problems.  In  general  if  the  Mach  number  is  high 
(>  4.0),  a  more  accurate  PEST  is  required.  Also, 
to  increase  the  stability  of  the  scheme  at  high 
Mach  numbers,  I0UT  should  be  read  in  as  1. 

Low  Mach  numbers  over  flared  bodies  also 
may  present  difficulties.  In  general,  the  Mach 
number  at  the  tip  of  the  base  should  not  be  less 
than  1.30.  In  addition,  special  care  should  be 
taken  such  that  the  Mach  number  along  the  flare 
does  not  become  subsonic. 

COMPUTER  PROGRAM  SUBROUTINES 

A  complete  list  of  the  subroutines  and  a  brief  description  of 
their  operation  are  given  below. 

KBPR.CNTROL 


KBPR.CNTROL  is  the  main  program  in  the  calculation  procedure 
(see  Flow  Chart  2  in  Appendix  D) .  Its  function  is  to  direct  the  other 
programs  in  the  correct  sequence.  Data  for  this  operation  of  the 
characteristics  solution  are  read  in  by  this  program.  In  addition, 
the  various  options  associated  with  the  characteristics  solution  are 
determined  here.  The  relative  simplicity  of  this  program  makes  it 
possible  to  perform  supplementary  calculations  and  include  additional 
steps  without  changing  the  more  complex  subprograms.  In  particular, 
this  program  directs  the  calculation  of  the  flow  field  over  the  body. 
The  direction  of  the  base  pressure  solution  is  controlled  by  the 
subroutine  KBPR.BASE5  which  is  called  by  KBPR.CNTROL. 

KBPR. FLARE 


The  subroutine  KBPR. FLARE  directs  the  actual  calculation  of  the 
characteristics  solution  both  over  the  body  (see  Figure  9)  and  in  the 
near  wake  region.  KBPR. FLARE  sets  up  the  grid  spacing  of  the  initial 
characteristic.  It  also  calculates  various  Mach  numbers  and  angles 
associated  with  the  base  pressure  solution. 
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KBPR.LCHAR 


The  subroutine  KBPR.LCHAR  determines  the  line  1  (or  2  )  as 

shown  in  Figures  4  and  8.  This  subroutine  is  also  responsible  for 
stopping  characteristic  calculations  beyond  this  line  so  that  useless 
computations  are  not  made.  The  variables  X2.  R2.  v2  and  02  are  the 
results  obtained  from  this  subroutine. 

KBPR.SIMR 

KBPR.SIMR  is  a  function  subprogram  which  performs  a  Simpson's  Rule 
integration  of  a  function  between  the  limits  of  XI  and  XN.  The 
integration  is  divided  into  N  equal  intervals  where  N  must  be  an  even 
number.  This  subprogram  was  developed  by  Professor  R.S.  Eikenberry  of 
the  University  of  Notre  Dame  Aero-Space  Department. 

KBPR.DINCR 

The  subroutine  KBPR.DINCR  performs  the  geometric  calculations 
which  determines  how  large  physically,  the  characteristics  solution  has 
to  be  such  that  there  are  sufficient  points  to  complete  the  calculations. 
This  returns  the  increment  of  spacing  between  grid  points  for  the  three 
locations:  along  the  body,  along  the  flare  or  boattail,  and  in  the  near 
wake  region.  These  spacings  are  functions  of  the  initial  Mach  number, 
the  flare  (or  boattail)  angle,  and  the  body  length. 

KBPR.AI1 


The  function  KBPR.AI1  is  the  argument  of  the  II  integral  used  in 
the  base  pressure  solution,  i.e,: 


4dn 


KBPR.AIl  = 


1  -  C2  <f>2 

3a 


KBPR.AI2 


KBPR.AI2  is  a  function  subprogram  which  defines  the  argument  of  the 
12  integral  used  in  the  base  pressure  solution,  i.e.: 

$2  dn 

KBPR.AI2  =  2  2 

3a 
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KBPR.AJ1 


y 


The  function  KBPR.AJ1  defines  the  argument  of  the  J1  integral  used 
in  the  base  pressure  solution,  i.e. : 

$  n  dn 

KBPR.AJ1  =  j  .  C2  ^2 


i 


KBPR.AJ2 

The  function  KBPRA.1"’  defines  the  argument  of  the  integral  J2 
which  is  used  in  the  base  pressure  solution,  i.e.: 


KBPR.AJ2 


‘i2  o  dn 


KBPR.I1 


The  functiov  KBPR. II  evaluates  the  II  integral  between  specific 
limits  usir,7  a  Simpson's  Rule  integration  (KBPR.SIMR). 

KBPR.I2 

KBPR.I2  is  a  function  subprogram  which  evaluates  the  12  integral 
using  Simpson's  Rule  integration. 

KBPR.J1 

KBPR.J1  is  a  function  subprogram  which  evaluates  the  J1  integral 
using  Simpson's  Rule  integration. 

KBPR.J2 

The  function  KBPR.J2  evaluates  the  J2  integral  using  Simpson's 
Rule  integration. 

KBPR . BODY 


The  subroutine  KBPR. BODY  calculates  the  distance  XBASE  (i.e.,  the 
axial  location  of  the  base) .  It  should  be  noted  that  this  calculation 
is  also  done  in  KBPR.CNTROL,  and  for  the  present,  KBPR. BODY  is  limited 
to  a  diagnostic  program;  it  tells  if  body  shapes  other  than  conical 
flares  or  boattails  are  used.  This  subprogram  is  still  included  since, 
in  the  future,  additional  afterbody  configurations  may  be  studied  which 
would  necessitate  such  a  subroutine. 
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KBPR.TAB 


This  function  subprogram  evaluates  a  single  function  for  a  given  , 

argument  by  interpolation  between  tabulated  values.  The  maximum  table  j 

size  is  16,384  entries,  and  the  maximum  degree  of  interpolation  is  ten.  f 

Table  search  is  by  a  binary  method  requiring  15  steps  to  find  any  entry 
in  a  table  of  16,384  entries.  Interpolation  (for  order  of  interpolation  f> 

greater  than  1)  is  done  by  means  of  Lagrange  Polynomials.  If  the  I 

argument  is  outside  the  limits  of  the  table,  the  first  (or  last) 
tabular  entry  is  returned  to  the  calling  program.  This  subprogram  was 
developed  by  Prof.  R.S.  Eikenberry  of  the  Aero-Space  Department  at  the 
University  of  Notre  Dame. 

KBPR.ERF 

The  function  subprogram  KBPR.ERF  evaluates  the  error  function  for 
a  given  argument. 

KBPR.PRES 

The  subroutine  KBPR.PRES  calculates  the  Mach  number,  pressure 
coefficient  (C^) ,  and  the  pressure  ratio  (P/PQ)  at  the  different 

X- locations  on  the  surface  of  the  body.  (It  should  be  noted  that  the 
characteristics  solution  assumes  an  inviscid  fluid  and,  therefore,  the 
no-slip  condition  on  the  body  surface  is  not  applicable.)  The  Mach 
number  is  obtained  by  the  use  of  the  subroutine  KBPR.PMTURN.  The 
pressure  ratio,  P/PQ,  is  determined  from: 

-y 

P/Po  ■  [‘  *  ¥  "2]  I 

The  pressure  coefficient  (based  on  the  free  stream  Mach  number)  is  i 

expressed  in  the  form: 

c  -  p/p°°  -  1  | 
p'iYM.!  ;| 

The  values  of  X,  M,  P/P  and  C  are  printed  out  at  each  location  in  : 

op 

KBPR.CNTROL. 
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KBPR.RSBF 


This  function  subprogram®  is  used  to  calculate  the  effective  sting 
radius  ratio  which  is  used  in  the  base  pressure  solution.  For  the 
present,  there  is  no  dependence  on  the  ratio  of  specific  heats,  y,  only 
the  Mach  number  at  the  tip  of  the  base.  In  the  near  wake,  the 
characteristics  solution  proceeds  until  a  grid  point  is  encountered 
which  is  less  than  this  sting  radius  (wake  radius).  In  KBPR.FBASE5,  a 
linear  interpolation  of  the  variables  is  made  to  this  sting  radius. 

KBPR. SIGMA 


This  function  subprogram8  calculates  the  jet  spread  parameter,  o. 
The  jet  spread  parameter  used  in  the  base  pressure  analysis  is  that 
formulated  by  Channapragada7 .  The  incompressible  value  of  the  jet 
spread  parameter,  oinc,  has  been  assumed  to  be  12.0. 


KBPR. LINEAR 

The  subroutine  KBPR. LINEAR  solves  by  successive  linear  interpolations 
the  value  of  XS  which  satisfies  FCN(XS)  =  0  +  ACX  where  ACX  is  an 
accuracy  requirement.  In  the  base  pressure  program,  FCN  =  KBPR.FBASE5 
=  C5  -  C6  and  XS  is  the  solution  for  P^/Pj.  XEST  is  the  initial 

estimate  of  the  base  pressure  ratio  (PEST),  and  DXLIN  is  the  size  of 
the  first  interval  of  the  independent  variable  (see  Figure  10). 


The  values  of  YEST  =  FCN(XEST)  and  Y1  =  FCN(XEST-DXLIN)  are 


:alculated.  The  slope,  S(= 


Y1  -  YEST 


),  is  used  to  project  toward  a  new 


DXLIN 

estimate  of  the  solution,  X2.  Then  Y2  =  FCN(X2)  is  evaluated.  If 
Y2  =  0  t  ACX,  then  X2  is  taken  as  the  solution  and  a  return  is  made  to 
the  calling  program.  If  not,  the  slope  is  then  taken  between  Y2  and  Y1 
and  this  new  slope  is  projected  again  toward  the  value  of  the  independent 
variable,  X3.  This  process  is  repeated  until  convergence  is  reached  by 
satisfying  the  accuracy  requirement. 


KBPR. CALC 


The  subroutine  KBPR. CALC  is  used  to  solve  for  the  characteristics 
solution  of  points  which  are  in  the  interior  flow  field  (i.e.,  not 
boundary  points).  The  procedure  for  calculating  the  location  of  a  new 


8.  Roaahe,  P.J.,  and  Mueller,  T.J.,  " A  FORTRAN  TV  Program  for 

Calculating  the  Axisyrmetric  Baee  Preeeure  in  Turbulent  Supersonic 
Flew  (BASED".  UNDAS  TN-866-M1,  University  of  Notre  Dame, 
Aerospace  Engineering,  August  1966. 
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grid  point  is  an  iterative  one  for  the  axisymmetric  method  of  character¬ 
istics  since  the  equation  is  no  longer  homogeneous.  Actually,  two 
iterations  are  used  and  this  yields  satisfactory  results  for  a  reasonable 
grid  spacing.  Once  the  X  and  R  coordinates  of  this  grid  point  have 
been  determined,  the  value  of  the  streamline  angle,  6,  and  the 
Prandtl-Meyer  turn  angle,  v,  are  calculated  for  that  point. 

The  two  upstream  points  which  are  used  to  calculate  this  new  point 
are  determined  by  the  subroutine  KBPR. FLARE  which  calls  for  KBPR.CALC. 

KBPR . SURF 


KBPR. SURF  is  the  subroutine  which  calculates  the  variables  X,  R, 
v  and  6  on  the  surface  of  the  body  (or  on  the  conetail  surface  in  the 
base  pressure  analysis) .  Again  this  becomes  an  iterative  procedure, 
and  two  guesses  are  made.  As  the  geometry  is  known,  it  is  only 
necessary  to  match  up  X  and  R  with  6,  as  6  at  the  surface  is  the  angle 
of  the  surface.  Once  this  has  been  done,  the  value  of  v  on  the  surface 
is  calculated.  Again  this  subroutine  is  directed  by  the  subroutine 
KBPR. FLARE. 

KBPR.CPB 


This  subroutine  provides  a  constant  pressure  boundary  base  pressure 
solution.  Again  this  is  an  iteration  procedure,  and  two  iterations  are 
made.  In  a  constant  pressure  boundary  solution,  the  value  of  v  remains 
unchanged.  Therefore,  the  correct  combination  of  X,  R  and  6  satisfy 
this  criterion.  This  subroutine  is  directed  again  by  KBPR. FLARE. 

It  should  be  noted  that  this  subroutine  is  not  called  for  a 
constant  pressure  boundary  solution  until  the  flow  has  passed  through 
the  Prandtl-Meyer  expansion  at  the  base.  To  calculate  the  characteristics 
through  this  expansion,  the  subroutine  KBPR. SURF  is  called  by  KBPR. FLARE. 

KBPR.SEPRTE 


This  subroutine  called  by  KBPR.CNTROL  determines  if  there  is  flow 
separation  on  the  surface  of  a  boattailed  afterbody.  The  separation 
criterion  which  is  used  compares  pressure  ratios  at  the  surface  with  a 
separation  pressure  ratio  for  that  free  stream  Mach  number  given  by 
Zukowski9  and  shown  in  Figure  11.  If  the  pressure  ratio  indicates 
separation,  a  message  is  printed  which  relates  the  pressure  ratio  which 
was  encountered  and  also  the  axial  distance  of  separation.  It  should  be 
noted  that  the  separation  criterion  used  is  valid  only  for  free  stream 
Mach  numbers  between  1.50  and  6.50.  For  Mach  numbers  not  within  this 
range,  the  separation  criterion  is  not  applied. 


9.  Zukowski,  £ . E ., "Turbulent  Boundary _Layer  Separation  in  Front  of  a 
Forward- Facing  Step".  AIAA  Journal,  Vol.  5,  No.  10,  pp.  1746-1753, 
October  1967, 
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KBPR . SETUP 


The  KBPR. SETUP  subroutine  is  used  in  an  attempt  to  locate  a  grid 
point  on  the  surface  of  the  body  which  is  as  close  to  the  end  of  the 
body  (or  region  of  the  body)  as  possible.  Referring  to  Figure  12a,  it 
is  noted  that  the  characteristics  solution  on  the  surface  has  not  yet 
exceeded  the  end  point,  XFIN.  However,  in  Figure  12b,  the  succeeding 
characteristic  solution  on  the  surface  has  exceeded  XFIN.  A  linear 
interpolation  is  now  made  between  X(I-2,1)  and  X(I,1)  such  that: 


R  TO 


X(I,1)  -  XFIN 
X(I , 1)  -  X(I2, 1) 


where  RTO  is  the  ratio  of  these  two  distances.  It  should  be  noted  that 
on  the  surface,  only  odd  values  of  I  appear  (see  Figure  6)  and, 
therefore,  X(I-2,1)  is  used  instead  of  X(I-1,1).  The  value  of  D1  which 
resulted  in  the  point  is  "lowered"  proportionally  to  RTO.  That  last 
characteristic  line  is  recalculated.  The  results  of  this  procedure  are 
shown  in  Figure  12c.  The  whole  procedure  is  again  repeated  to  achieve 
further  accuracy.  The  KBPR. SETUP  subroutine  then  makes  normal  return 
to  the  subroutine  KBPR. FLARE  and  the  remainder  of  the  characteristic 
lines  are  calculated. 


KBPR. PMTURN 


This  subroutine  solves  the  inverse  Prandtl-Meyer  turn  equation  of 
an  ideal-gas  compressible  flow.  It  either  returns  the  downstream  Mach 
number  after  an  isentropic  turn,  or  simply  determines  the  Mach  number 
for  a  given  Prandtl-Meyer  turn  angle,  v,  where  v  is  defined  as: 

v  =  J\ ~T  tan":  J\ T'T  CM'2'  ~  1)  -  fan’1  M7  -'  1 

Calculations  are  made  using  a  Newton-Raphson  technique. 

KBPR.PMANGL 

KBPR.PMANGL  is  a  function  subprogram  which  solves  the  Prandtl-Meyer 
turn  equation  of  an  ideal-gas,  compressible  flow.  It  returns  the 
Prandtl-Meyer  turn  angle,  v,  which  is  defined  by  the  Equation  above. 

KBPR.S0TE2B 


The  subroutine  KBPR.S0TE2B  solves  a  double-valued  trancendental 
equation  FCN(XS)  =  0  ±  AC,  where  AC  is  an  accuracy  requirement.  In  the 
base  pressure  analysis,  FCN  is  the  function  FSTEP4(ETAJ3) .  This  is  a 
complicated  integral  equation  to  be  solved  for  n^3  =  ETAJ3,  the  upper 

limit  of  integrals  12  and  J2.  To  use  S0TE2B,  the  upper  extreme,  XU, 
and  the  lower  extreme,  XL,  must  be  defined.  The  values  of  XU  and  XL 
have  been  set  at  +3  and  -3,  respectively.  The  indicator,  IND  =  +1,  causes 
the  upper  solution  to  be  given,  while  IND  <  0  gives  the  lower  solution. 
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Figure  13a  shows  the  qualitative  behavior  of  FSTEP4  and  illustrates 
the  method  of  solution.  The  range  of  X  for  a  solution  is  first  reduced 
to  a  AX  which  contains  the  upper  solution  XS  (for  IND  =  +1),  but  not 
the  lower  solution.  To  do  this  a  search  is  initiated  for  a  sign  change 
in  FSTEP4  as  X  is  decreased  from  XU  =  +3.  To  avoid  possibly  skipping 
over  the  sign  changes,  the  technique  shown  in  Figure  13b  is  used.  The 
sweeps  continue  until  a  sign  change  is  found  between  some  XI  and  X2. 

Now  the  problem  is  reduced  to  finding  a  solution  XS  to  a  single-valued 
trancendental  equation  in  which  XS  is  known  to  lie  between  XI  and  X2. 

The  solution  is  found  by  interval  halving  which  is  shown  in  Figure  13c. 
The  calculation  proceeds  A-B-C-D-E  in  which  the  interval  is  decreased 
by  half  at  each  trial.  This  technique  is  vejjy  rapid  as  n-calculations 
reduce  the  interval  sice  to  less  than  (1/2)  of  its  initial  size. 

KBPR. SHOCK 


The  subroutine  KBPR. SHOCK  gives  a  closed-form  solution  of  the 
shock  wave  cubic  equation  for  an  ideal  gas.  It  returns  to  the  calling 
program  the  wave  angle  when  the  initial  Mach  number  and  streamline  turn 
angle  are  given.  If  the  detachment  turn  angle  is  exceeded,  the  value 
returned  is  the  detachment  turn  angle.  The  indicator,  N=4,  provides  a 
weak  shock  solution,  while  S-0  gives  the  strong  shock  solution.  The 
form  of  the  exact  solutions  may  be  found  in  Reference  7.  This  program 
was  developed  by  Patrick  J.  Roache  of  the  Aero-Space  Department  at  the 
University  of  Notre  Dame. 

KBPR. FSTEP4 

The  function  subprogram  KEPR.FSTEP4  is  the  subprogram  in  the  base 
pressure  section  of  the  program  which  combines  the  integrals  resulting 
in  the  calculation  of  the  geometric  parameter.  The  difference  between 
the  actual  geometric  parameter  and  the  calculated  geometric  parameter 
is  what  is  transferred  into  KBPR.S0TE2B. 

KBPR . BASE5 


The  subroutine  KBPR.BASE5  is  the  program  which  directs  the 
calculation  of  the  base  pressure  at  the  rear  of  the  body.  This 
subroutine  is  called  by  KBPR.CNTROL  after  the  flow  field  over  the  body 
has  been  established.  In  KBPR. BASES,  the  parameters  for  the  operation 
of  the  base  pressure  section  of  the  program  are  read  in.  This  subroutine 
is  primarily  used  for  reading  data  and  printing  out  the  variables  once 
a  solution  has  been  obtained. 

KBPR. FBASE5 


This  function  subprogram  is  called  by  KBPR.BASE5  and  is  the 
directing  program  in  the  base  pressure  analysis.  The  effective  sting 
radius  (see  KBPR.RSBF)  is  determined  and  then  the  geometric  parameter 


24 


DISCUSSION  OF  RESULTS  FOR  CONICAL  AFTERBODIES 


A  large  amount  of  experimental  afterbody  and  base  pressure  data 
have  been  taken  for  boattailed  bodies  of  revolution.  These  investi¬ 
gations  include  conical,  circular  arc  and  parabolic  shaped  boattails 
of  various  length  to  diameter  ratios  for  Mach  numbers  from  about  1.7 
to  about  6.0.  Because  of  the  possible  occurrence  of  flow  separation 
as  well  as  model  support  or  sting  interference,  these  data  must  be 
examined  very  carefully  before  comparisons  can  be  made  with  the 
analytical  results. 

The  computer  program  BOATTAIL/FLARE-BASE5  can  be  used  to  calculate 
the  turbulent  base  pressure  for  configurations  similar  to  all  those 
shown  in  Figure  1.  Both  non-uniform  and  uniform  flows  approaching 
the  base  or  afterbody  can  be  handled.  Sample  inputs  with  the 
corresponding  outputs  are  included  in  Appendix  C. 


UNIFORM  FLOW  APPROACHING  THE  AFTERBODY 

This  computer  program  has  been  used  and  the  results  verified  for 
uniform  approach  Mach  numbers  from  1.3  to  5.0,  ratios  of  specific  heats 
from  1.2  to  1.667,  base  temperature  ratios  from  0.5  to  3.0,  and 
conical  boattail  and  flare  angles  up  to  15°.  A  few  examples  of  the 
type  of  results  obtained  fr*.  this  program  will  be  presented.  All 
analytical  results  in  this  report  are  for  the  case  of  rising  pressure 
jet  mixing. 

The  results  obtained  from  the  B0ATTAIL/FLARE-BASE5  computer 
program  for  y  =  1.4,  T./T  =  1.0  and  an  initial  uniform  Mach  number  of 

2.0  are  shown  in  Figure  14.  The  geometry  of  the  boattail  body  studied 
analytically  was  the  same  as  that  studied  experimentally  by  Reid  and 
Hastings10  so  that  a  direct  comparison  could  be  made.  The  experiments 
of  Reid  and  Hastings  closely  approximate  the  analytical  flow  model  used 
i.e.,  uniform  external  flow  except  for  a  relatively  thin  turbulent 
boundary  layer.  Other  data  available  either  include  nose  effects,  angle 
of  attack  of  effects  and/or  an  undetermined  boundary  layer  character  and 
thickness  approaching  the  base.  Excellent  agreement  between  the  present 
analytical  method  and  the  data  of  Reid  and  Hastings  is  indicated  in 
Figure  14  for  r^/r  =  0.58.  The  computer  solution  was  unable  to 


10.  Reid ,  J. ,  and  Hastings,  R.C.,  "Experiments  on  the  Axi-Symmetric 
Flaw  Over  Afterbodies  and  Bases  at  M  =  2. 0".  Royal  Aircraft 
Establishment  Report  No.  Aero. 2628,  October  1959, 

(STAR  NB6- 20065). 


produce  reliable  results  for  boattail  angles  (-4)  lower  than  about  3° 
because  of  the  grid  spacing  of  the  characteristics  solution  necessary 
to  achieve  accuracy.  For  small  angles  and  moderate  r^/r^.,  the  boattail 

(or  flare)  length  becomes  quite  long,  requiring  a  large  grid  size. 
However,  if  the  curve  of  Figure  14  is  extrapolated  to  a  boattail  angle 
of  zero,  the  resulting  base  pressure  ratio  agrees  with  the  results  for 
the  flow  over  a  cylindrical  body  obtained  by  Mueller3,  i.e.,  P^/Pj  =  0.6. 

Results  for  a  7°  boattail  and  y  =  1.4,  T^/T^  =  1.0,  L/rc  =  0.75  for 

Mach  numbers  from  1.3  to  5.0  are  shown  in  Figure  15. 

The  influence  of  afterbody  angle,  4,  on  the  base  pressure  ratio 
for  initial  uniform  Mach  numbers  2.0,  3.0  and  4.0  is  presented  in 
Figure  16.  These  results  were  obtained  for  a  fixed  afterbody  length, 

L/rc  =0.5,  a  uniform  initial  Mach  number,  y  =  1.4,  and  T^/T^  =  1.0. 

Figure  17  shows  the  influence  of  initial  Mach  number  on  the  base 
pressure  ratio  for  several  afterbody  angles  and  a  longer  afterbody.  In 
this  figure  a  fixed  afterbody  length  of  L/r^  =  2.0,  i.e.,  a  one  caliber 

afterbody,  was  used  with  a  uniform  initial  Mach  number,  y  =  1.4,  and 

T./T  =  1.0.  These  results  indicate  the  well  behaved  nature  of  the 
b  oa 

analytical  results  for  conical  afterbodies.  If  desired,  this  method 
could  be  used  to  study  the  influence  of  afterbody  shape,  y,  and  T./T 

u  03. 

on  the  base  pressure  ratio.  A  sample  input  and  output  from  BOATTAIL/ 
FLARE-BASES  for  uniform  flow  of  M  =  3.0  approaching  a  one  caliber  7° 
boattail  is  shown  at  the  end  of  Appendix  C. 


NON-UNIFORM  FLOW  APPROACHING  THE  AFTERBODY 

A  large  amount  of  effort  has  been  expended  to  study  the  use  of 
B0ATTAIL/FLARE-BASE5  for  non-uniform  supersonic  flows  approaching  the 
afterbody.  For  one  of  these  studies  a  linear  Mach  number  distribution 
approaching  a  conical  boattail  with  L/rc  =0.5  was  chosen.  The 

non-uniform  Mach  number  profile  was  input  at  the  end  of  the  cylindrical 
body  and  ranged  from  M£  =  2.10  at  the  body  surface,  i.e.,  r/rc  =  1.0, 

to  M  =  2.3  at  r/r  =  1.6  from  the  body  surface.  The  freestream  at 
oo  c 

M^  =  2.3  was  continued  to  r/rc  =  2.0.  The  conical  boattail  angle  was 

varied  from  0°  to  9°.  A  comparison  of  the  base  pressure  ratio  calculated 

for  a  uniform  flow  with  M  =2.3  and  the  linear  non-uniform  flow  is 

c 

presented  in  Figure  18  as  a  function  of  boattail  angle.  This  non-uniform 
Mach  number  profile  produces  base  pressure  ratios  which  are  lower  than 
the  values  obtained  for  uniform  flow. 
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To  compare  the  base  pressure  ratios  calculated  using  BOATTAIL/ 
FLARE-BASE5  with  experimental  values  for  a  particular  U.S.  Army 
configuration,  several  specific  cases  were  run.  The  configuration 
studied  had  a  three  caliber  secant  ogive  nose,  a  two  caliber  cylindrical 
section  and  a  one  caliber  7°  conical  boattail.  The  non-uniform  Mach 
number  distribution  just  upstream  of  the  boattail,  obtained  from  a  BRL 
solution  of  the  Euler  Equations  using  MacCormack's  predictor- corrector 
numerical  technique1!.  This  non-uniform  input,  along  with  the  correspond¬ 
ing  output,  is  shown  at  the  end  of  Appendix  C.  Figure  19  presents  a 
comparison  of  the  computed  results  for  this  non-uniform  flow  input  as 
well  as  for  many  uniform  flow  cases  with  the  BRL  experimental  data. 

As  the  non-uniform  Mach  number  distribution  was  close  to  a  uniform 
profile,  it  is  not  surprising  that  this  result  falls  on  the  uniform 
flow  curve.  The  fact  that  the  experimental  data  for  the  highest 
Reynolds  number  tested  also  fall  on  this  curve  is  an  indication  of  how 
well  base  pressure  ratio  is  predicted  by  B0ATTAIL/FLARE-BASE5 . 

It  was  found  by  L.D.  Kayser  of  BRL  that  this  program  could  also 
be  used  to  obtain  base  pressure  results  for  cones.  By  choosing  a 
fairly  long  flare  with  an  angle  equal  to  the  cone  angle  and  a  uniform 
Mach  number,  results  which  agreed  well  with  experimental  data  were 
obtained. 


CONCLUSION  AND  RECOMMENDATIONS 

The  computer  program  now  operational  on  the  Aberdeen  Proving 
Ground  UNIVAC  1108  has  been  found  to  produce  results  which  agree  well 
with  experimental  data  for  high  Reynolds  numbers.  These  programs  will 
run  for  a  wide  variety  of  geometries  and  approach  flow  conditions  and 
can  even  be  used  to  calculate  base  pressure  on  cones.  These  programs 
can  be  modified,  if  desired,  to  include  the  effects  of  base  bleed  and 
initial  boundary  layer  thickness. 


11.  Sander's,  B.R. ,  "Three-Dimensional,  Steady  Invisoid  Flow  Field 
Calculations  with  Application  to  the  Magnus  Problem",  PhD 
Dissertation,  University  of  California,  Davie,  California,  May  1974. 


Figure  2.  Flow  Model  for  the  Influence  of  Initial  Flow  Direction  on  the 
Axisymmetric  Turbulent  Base  Pressure. 


Figure  3.  Trailing  Wake  Radius  Ratio  Versus 
Mach  Number  for  y  =  1.4  (Ref.  6) 
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Figure  4.  Schematic  of  Axi symmetric  Characteristic  Grid 


Figure  5b.  Characteristic  Construction  Technique  (I  =  Even  Number) 


Characteristic  Numbering  Scheme 


igure  7.  Characteristic  Termination  Lines 


Figure 


Variation  of  X2,  R2,  v2,  e2  along  a  Characteristic  Line, 
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Interpolation  Scheme  Used  in  Subroutine  LINEAR 


Figure  11.  Variation  of  Separation  Pressure  Ratio 
with  Mach  Number  from  Zukowski  (Ref.  9) 
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XFIN 


Figure  12a.  Characteristics  Mesh  Before  Reaching  XFIN 


Figure  12b.  Characteristics  Mesh  After  Reaching  XFIN 


Figure  12c.  Corrected  Characteristics  Mesh  Done  by  Subroutine  SETUP 


Figure  13a.  Qualitative  Behavior  of  Function  Subprogram  FSTEP4 


XXX  STEP 

X  X  STEP  '2 

X  X  X  X  STEP  -3 

XXXXXXXX  STEP  *4 
XXXXXXXXXXXXXXXX  STEP  -5 

Figure  13b.  Searching  Scheme  Used  to  Locate  the  Regions 
of  Solution  in  S0TE2B 


Figure  15.  Base 


I 


Pressure  Ratio  Versus  Boattail  Angle  for  Uniform 


the  Calculated  Base  Pressure  Ratio  for  Uni 
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LIST  OF  SYMBOLS 


Variable  defined  by  Eqn.  2 


Crocco  number 


Pressure  coefficient  base  on  M 

00 

Distance  along  last  characteristic  line 

Distance  along  initial  characteristic  line 

Subscript 

Subscript 

Afterbody  length 

Total  body  length 

Mach  number 

Direction  perpendicular  to  streamline 
Absolute  pressure 

Reference  streamline  close  to  but  outside  of  the 
mixing  region;  radial  location  of  characteristic 
grid  point;  empirical  compressible  divergence  factor 

Radial  location  of  grid  point  along  last  characteristic 

Radius 


Radius 


Direction  along  a  streamline 
Temperature 

Velocity  in  x  of  X  direction 
Velocity  along  streamline 


Cordinates  of  the  reference  (inviscid)  coordinate 
system 
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X  Axial  location  of  grid  points  in  characteristic  grid 

x,  y  Coordinates  of  the  intrinsic  (viscous)  coordinate 

system 

X2  Axial  location  of  grid  point  along  last 

characteristic  line 

a  Dummy  variable  in  Eqn.  14 

8  Temperature  ratio  T0a/Tb 

y  Ratio  of  specific  heats 

5  Distance  along  right  running  wave 

n  Distance  along  left  running  wave;  similarity 

parameter 

0  Streamline  angle 

02  Streamline  angle  at  grid  point  along  last 

characteristic  line 

u  Mach  angle 

v  Prandtl-Meyer  turn  angle 

v2  Prandtl-Meyer  turn  angle  at  grid  point  along  last 

characteristic  line 

a  Jet  spread  parameter 

<t  Afterbody  angle 

<j>  Dimensionless  velocity  defined  as  u/ua 


\ 

I 


SUBSCRIPTS 

1,2, 3, 4  Refer  to  conditions  at  cross-sections  of  basic 
flow  model 

a  Refers  to  conditions  in  the  external  stream 

adjacent  to  the  mixing  region 

b  Refers  to  condition  at  the  base 
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c  Refers  to  the  conditions  upstream  of  flare  or 

boattail 

d  Refers  to  conditions  along  a  streamline  whose  kinetic 

energy  is  just  sufficient  to  enter  the  recompression 
region 

inc  Incompressible  value 

j  Condition  along  the  jet  boundary  separating 

streamline 

m  Coordinate  shift  in  the  mixing  theory  due  to  the 

momentum  integral 

o  Refers  to  stagnation  conditions 

R  Condition  along  the  R  streamline 

s  Refers  to  separation  conditions 

w  Refers  to  wake 

“  Refers  to  freestream  conditions 

(  )  Refers  to  iteration  number 

SUPERSCRIPTS 

—  Average  value 


APPENDIX  A 


SYMBOL  DEFINITIONS 

The  symbols  used  in  each  subroutine  are  defined  under  each 
subroutine  title.  The  flow  model  used  in  the  axisymmetric  base 
pressure  analysis  (Fig.  2)  and  the  afterbody  configurations  (i.e. 
boattail  or  flare  shown  in  Fig.  9)  are  necessary  references  for 
many  definitions. 

KBPR.CNTROL 

BS  =  Afterbody  shape 

Cl  =  [X2(I)-X2(I-1) r 

C2  =  [ R2 ( I) -R2 (1-1) ] 2 

CDEL1B  =  Determines  distance  from  body  of  first  point  in 
characteristic  field  (used  at  the  base) 

CDEL1P  =  Determines  distance  from  body  of  first  point  in 
characteristic  field  (used  along  the  body) 

CP  =  Pressure  coefficient 

D  =  Distance  along  a  characteristic  (=  /Cl  +  C2  ) 

DEL 12  =  XFLARE  -  XSTART 

DEL23  =  XBASE  -  XFLARE 

DINC1  =  Used  to  calculate  maximum  spacing  of  grid  points  in  the 
characteristic  field  for  the  cylindrical  surface 
DINC2  =  Used  to  calculate  maximum  spacing  of  grid  points  for 
afterbody  region 

DINCB  =  Used  to  calculate  maximum  spacing  of  grid  points  for  the 
base  region 

G  =  Ratio  of  specific  heats  (=  Cp/Cv) 

Gi  =  -G/(G-1) 

I  =  Subscript 

IBOUND  =  Determines  type  of  base  pressure  solution  desired 
ILIMA  =  Maximum  number  of  grid  points  in  X-direction 
IOPTL  =  Option  number  permitting  variation  of  input  data 


IOUT 


=  Determines  how  characteristics  solution  will  be  trans¬ 
ferred  to  base  region 
IREAD1  =  Determines  how  initial  profile  will  be  set  up 

IREAD2  =  Determines  how  characteristics  solution  will  be  trans¬ 
ferred  to  the  afterbody  region 
JL1MA  *  Maximum  number  of  grid  points  in  R-direction 

JLIMT  *  Number  of  points  to  be  read  in  for  a  non-uniform  profile 

M  «  Mach  number  in  characteristics  field 

MACH  =  Mach  number  along  initial  profile 

MESHPM  =  Number  of  discrete  turns  into  which  a  Prandtl-Meyer 

expansion  is  divided 

MIA  ■  Mach  number  at  beginning  of  characteristics  solution  at 
each  section 

NU  *  Prandtl-Meyer  turn  angle  in  characteristics  field 

NU2  *  Prandtl-Meyer  turn  cngle  along  initial  profile 

OMACH  *  Reference  Mach  number  (usually  M  ) 

00 

PBPINF  «  Pb/P,„ 

PBPO  "  Pb/P0 

PBP1  =  Pb/Pi 

PHI  =  Flare  (or  boattail)  angle  (radians) 

PHID  ■  PHI  (degrees) 

P1NF  *  P„/P0 

PSURF  =  P/PQ  al°n8  the  surface  of  the  body 

P10P0  =  Px/P0 

R  =  Radius  in  characteristic  field 

RBASE  ■  Base  radius  (non-dimensionalized) 

RBODY  »  Body  radius  (non-dimensionalized) 

RBODYA  »  Actual  body  radius 
RBSRSF  «  RBASE /RBODY  or  XBASE 

RSRB  »  Rs/Rbase  where  Rs  *  theoretical  sting  radius 

R2  ■  Radius  along  initial  profile 

SMACH  »  Surface  Mach  number  along  body 
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T  =  Streamline  angle  (radians)  in  characteristics  field 

TS  =  T( 1 , 1) 

T2  =  Streamline  angle  of  initial  profile 

X  =  Axial  distance  in  characteristics  field 

XBASE  =  Distance  of  base  from  origin 

XFLARE  =  Distance  to  the  beginning  of  the  afterbody 

XSTART  =  Distance  to  the  beginning  of  flow  field  calculations 

XSURF  =  Distance  along  the  body  surface 

X2  ■  Axial  distance  of  initial  profile 

KBPR. FLARE 

CAXGLE  =  EHUIA  -  T(l,l) 

CDEL1  =  Determines  distance  from  body  of  first  point  in  charac¬ 
teristic  field 
Cl  =  [X2(I)-X2(1-1) ]2 

C2  =  [R2(I)-R2(I-1) ]2 

C3AS  *  (C3a) 2 

D  =  Distance  along  a  characteristic  (=•  /Cl  +  C2  ) 

DEGNU  =  Prandtl-Meyer  turn  angle  (degrees) 

DEGT  =  Streamline  angle  (degrees) 

DEL  =  Increment  size  along  initial  characteristic 

DELMX  =  Maximum  increment  size  along  initial  characteristic 

DELI  =  Size  of  first  increment  along  initial  characteristic 

DINC  =  Used  in  the  calculation  of  DELMX 

DTC  =  TS/NTC 

D1  *  Distance  along  initial  profile 

EJK  =  Used  to  increment  J 

KKJ  =  Used  to  increment  J 

EMFSRB  *  Mach  number  at  position  3a  in  base  pressure  model 

EMU1A  »  Mach  number  at  tip  of  base 

EM2A  =  Mach  number  at  position  2a  in  base  pressure  model 

1  »  Subscript  in  axial  direction 

IA  =  Option  number  used  in  subroutine  KBPR. SETUP 
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I BOUND 
IC 

ICORN 

ICRNPT 

II 

III 

ILIMA 

J 

JCORN 

JLIMA 

JLIMT 

M 

MESHPM 

NTC 

NTCMIN 

NU 

NU2 

R 

RBODYA 

RLAST 

RMOD 

RSRB 

R2 

R2P 

SINANG 

SOL3A 

T 


Determines  type  of  base  pressure  solution  desired 
1-1 

Number  of  values  of  I  necessary  to  turn  through  an 
expansion  at  the  surface 
ICORN  +  20 
1/2 

Determines  body  location  (*  1  along  cylinder  part  of 
body;  =  2  along  afterbody;  *  3  in  base  flow  region) 
Maximum  number  of  grid  points  in  X-direction 
Subscript  in  radial  direction 

Number  of  values  of  J  necessary  to  turn  through  an 
expansion  at  the  surface 

Maximum  number  of  grid  points  in  R-direction 
Number  of  points  to  be  read  in  for  a  non-uniform  profile 
Mach  number  in  characteristics  field 
Number  of  discrete  turns  into  which  a  Prandtl-Meyer 
expansion  is  divided 

Number  of  turns  into  which  a  Prandtl-Meyer  expansion  is 
divided 

Minimum  value  of  NTC 

Prandtl-Meyer  turn  angle  in  characteristics  field 
Prandtl-Meyer  turn  angle  along  initial  profile 
Radius  in  characteristic  field 
Actual  body  radius 
R(L, 1) 

Diraentionalized  radius 

Rg/Rbase  where  Rs  *  theoretical  sting  radius 
Radius  along  initial  profile 

RCl.l) 

Sin  (CANGLE) 

M3 A/ M2 A  (or)  THET12 

Streamline  angle  (radians)  in  characteristics  field 
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THET12 

THET3A 

TS 

T2 

V3A 

X 

XMOD 

X2 

X2P 


Cl 

C2 

D 

DEGNU 

DEGT 

EKJ 

III 

II 

13 

JLIMT 

J1 

M 

NU 

NU2 

R 

RATIO 

RBODYA 

RMOD 

R2 

T 


Angle  of  streamline  off  base 
Angle  of  streamline  at  recompression 
T(1 ,1) 

Streamline  angle  of  initial  profile 

Prandtl-Meyer  turn  at  Rs  (Rs/Rbase  where  Rs  *  theoretical 
sting  radius) 

Axial  distance  in  characteristics  field 
Dimensionalized  axial  distance 
Axial  distance  of  initial  profile 

X(l,l) 

KBPR.LCHAR 

[X2(I)-X2(I-1)]2 

[R2(I)-R2(I-1)32 

Distance  along  a  characteristic  (■  /Cl  +  C2  ) 
Prandtl-Meyer  turn  angle  (degrees) 

Streamline  angle  (degrees) 

Used  to  increment  J 

Location  on  body  (=  1  along  cylindrical  part  of  body; 

*  2  along  afterbody;  =  3  in  base  flow  region) 

Value  of  I  at  last  point  in  region  1 
I  +  II-l 

Number  of  points  to  be  read  in  for  a  non-uniform  profile 
1-1 

Mach  number 

Prandtl-Meyer  turn  angle  in  characteristics  field 
Prandtl-Meyer  turn  angle  along  initial  profile 
Radius  in  characteristic  field 

Used  to  determine  closest  point  to  the  end  of  a  region 
Actual  body  radius 
Dimensionalized  radius 
Radius  along  initial  profile 

Streamline  angle  (radians)  in  characteristics  field 
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T2 

X 

XMOD 

X2 

DX 

I 

N 

SEV 

SOD 

SIMR 

X 

XI 
XN 


«  Streamline  angle  of  initial  profile 
=  Axial  distance  in  characteristics  field 

*  Dimensionalized  axial  distance 

=  Axial  distance  of  initial  profile 

KBPT.SIMR 

=  An  increment  along  the  abscissa 
=  Subscript 

*  Number  of  increments  (even  number) 

=  Sum  of  the  even  terms 

=  Sum  of  the  odd  terms 
=  Simpsons  Rule  integration  result 
=  Value  of  independent  variable 
=  Lower  limit  of  X 
=  Upper  limit  of  X 


ANG1 

ANG2 

ANG3 

DIFF 

DINC1 

DINC2 

DXNCB 

HI 

H1A 

H1EFF 

H2 


KBPT.DINCR 

*  MU-T1 

=  MU-T2 

=  MU-T3 

=  JLIMA-23 

=  Used  to  calculate  maximum  spacing  of  grid  points  in  the 
characteristic  field  for  the  cylindrical  surface 

*  Used  to  calculate  maximum  spacing  of  grid  points  for 
afterbody  region 

=  Used  to  calculate  maximum  spacing  of  grid  points  for  the 
base  region 

*  Additional  height  over  regions  2  and  3  which  is  needed 
such  that  a  characteristics  solution  will  be  complete 

*  Total  height  required  at  XSTART 

»  Effective  height  at  XSTART 

«  Additional  height  over  region  3  which  is  needed  such 
that  a  characteristic  solution  will  be  complete 
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H2A  =  Total  height  required  at  XFLARE 
H2EFF  »  Effective  height  at  XFLARE 

H3  *  Height  which  is  needed  such  that  a  characteristics  solu¬ 

tion  will  be  complete  in  the  base  region 
H3A  *  H3  +  RSTING 

H3EFF  *  Effective  height  at  XBASE 
LI  =  XFLARE-XSTART 

L2  =  XBAS  E-XFLARE 

L3  =  Length  from  the  base  to  the  sting  radius  along  a  5°  angle 

MU  =  Mach  angle 

RBASE  »  Base  radius 

RSTING  =  Theoretical  sting  radius 

T1  =  A  certain  percentage  of  the  angle  PHI  depending  on 

whether  PHI  is  positive  or  negative  (Region  1) 

T12  =  Lower  limit  of  flow  separation  angle  from  base  («  5°) 

T2  =  A  certain  percentage  of  the  angle  PHI  depending  on 

whether  PHI  is  positive  or  negative  (Region  2) 

T3  •  A  certain  percentage  of  the  angle  PHI  depending  on 

whether  PHI  is  positive  or  negative  (Region  3) 

W  *=  Used  to  convert  degrees  to  radians 

XBASE  «  Axial  distance  to  base  (see  Fig.  9) 

KBPR.AI1 

A  =*  Velocity  profile  [■*  -|  (1  +  erf  x)  ] 

All  -  Argument  of  II  integral 

2 

CS  »  (Crocco  number) 

TRB  «  Base  temperature  ratio  (-  Tb/T0a) 

TRF  -  A  +  TRB  (1-A) 

KBPR.AI2 

A  »  Velocity  profile  [“  y  (1  +  erf  x)] 

AI2  »  Argument  of  12  integral 

AS  -  A2 


cs 

= 

2 

(Crocco  number) 

TRB 

= 

Base  temperature  ratio  (=  Tb/Toa) 

TRF 

A  +  TRB  ( 1-A) 

KBPR.AJ1 

A 

= 

Velocity  profile  [«  —  (1  +  erf  x) ] 

AJI 

- 

Argument  of  J1  integral 

CS 

= 

2 

(Crocco  number) 

TRB 

= 

Base  temperature  ratio  (=  T[j/Toa) 

TRF 

A  +  TRB  (1-A) 

KBPR.AJ2 

A 

= 

Velocity  profile  [=  -|  (1  +  erf  x)  ] 

AJ2 

= 

Argument  of  J2  integral 

AS 

= 

A2 

CS 

= 

2 

(Crocco  number) 

TRB 

Base  temperature  ratio  (=  Tj,/Toa) 

TRF 

A  +  TRB  (1-A) 

KBPR.I1 

11 

Integral  11 

KBPR.I2 

12 

Integral  12 

KBPR.J1 

J1 

Integral  J1 

KBPR.J2 

J2 

Integral  J2 

KBPR . BODY 

DELTAR 

= 

^base  "  ^body 

DELTAX 

3 

^base  “  ^flare 

KBPR. TAB 

D1 

« 

Increment  of  I 

60 


DIL  =  Determines  how  many  points  below  I  will  be  needed 

DIU  »  Determines  how  many  points  above  I  will  be  required 

I  *  Subscript 

IL  -  I  -  DIL 

IU  =  I  +  DIU 

J  =  Subscript 

P  =  Lagrange  coefficient 

KBPR.ERF 

A1,A2 ,A3,A4,A5  -  Coefficients 

DENOM  =  Denominator  of  term 

ERF  =  Error  function 

y  =  Point  at  which  erf  is  to  be  evaluated 

V  =  |x| 

KB PR. PRES 

CP  =  Pressure  coefficient 

G1  =  -G/G-l) 

G2  -  (G-l)/2 

G  *  Ratio  of  specific  heats 

I  *  Subscript 

ILAST  «  I 

ILAST  1  =  Value  of  I  at  last  point  on  the  surface  in  Region  1 

ILAST  2  =  Value  of  I  at  last  point  on  the  surface  in  Region  2 

II  -  I 

12  =  1/2  +  1 

14,16  =  Dummy  integer  numbers 

M  «  Mach  number 

NU  =  Prandtl-Meyer  turn  angle 

PINFIN  -  Po/P*, 

PSURF  =  P/P0  along  the  surface  of  the  body 

R  =  Radius  in  characteristic  field 

RE1,  RE2,  RE3,  RE4 ,  RES,  RE6  ”  Dummy  real  numbers 


RTO  =  Used  to  determine  last  point  on  surface  at  each  region 
SM  =  Surface  Mach  number 

SMACH  =  Surface  Mach  number 

T  =  Streamline  angle 

X  =  Axial  distance 

XSURF  =  Axial  distance  of  surface 

KBPR.RSBF 

Ml  =  Mach  number  at  base 

RSBF  =  Function  which  describes  Rsting/^base 

KBPR. SIGMA 

2 

CS  =  (Croc co  number) 

R  =  Constant  (or)  function  of  CS 

SIGMA  =  let  spread  parameter 

SR  =  Jet  spread  parameter  ratio  (=  SIGMA/SIGMAj^qOMP^ 

KBPR. LINEAR 

ACX  =  Accuracy  desired 

DX  =  Increment  of  X 

S  =  Slope 

XI  =  ValueofXatpointl 

X2  =  Value  of  X  at  point  2 

X3  *  Projected  value  of  X 

XS  =  Solution  of  equations 

Y1  =  Value  of  Y  at  point  1 

Y2  =  Value  of  Y  at  point  2 

Y3  =  Value  of  Y  at  point  3 

KBPR . CALC 

Cl  =  Inhomogeneous  part  of  characteristic  equation  at  point  1 

C13  =  Projected  value  of  C  from  1  to  3 

C2  =  Same  as  Cl  except  at  point  2 

C23  =  Same  as  C13  except  from  2  to  3 

C3  =  Final  value  of  inhomogeneous  part  of  equation 
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ETA23 

- 

Distance  along  left  running  wave 

FI 

= 

Nul  +  T1  +  C13  x  Z13 

FIG 

= 

First  guess  at  FI 

F2 

= 

Nu2  -  T2  +  C23  x  ETA23 

F2G 

= 

First  guess  at  F2 

MUl 

= 

Mach  angle  at  1 

MU2 

= 

Mach  angle  at  2 

MU3G 

= 

First  guess  at  Mach  angle  at  3 

Ml 

= 

Mach  number  at  1 

M2 

= 

Mach  number  at  2 

M3G 

= 

First  guess  at  Mach  number  at  3 

NU3G 

= 

First  guess  at  Prandtl-Meyer  turn  angl 

R13 

= 

Average  radius  between  1  and  3 

R2  3 

= 

Average  radius  between  2  and  3 

SI 

= 

Slope  of  right  running  wave  at  1 

S2 

S 

Slope  of  left  running  wave  at  2 

T3G 

S 

First  guess  at  the  streamline  angle  at 

Z13 

Distance  along  right  running  wave 

KBPR. SURF 

Cl 

= 

(See  KBPR.CALC) 

Cl  3 

= 

(See  KBPR.CALC) 

DTC 

= 

(See  KBPR. FLARE) 

I 

= 

Subscript 

II 

- 

1/2 

MUl 

= 

(See  KBPR.CALC) 

MU3G 

- 

(See  KBPR.CALC) 

Ml 

B 

(See  KBPR.CALC) 

M3G 

B 

(See  KBPR.CALC) 

NTC 

S 

(See  KBPR. FLARE) 

NU3G 

- 

(See  KBPR.CALC) 

R13 

- 

(See  KBPR.CALC) 

R2 

- 

Radius  at  point  2 
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SI 

=s 

(See  KBPR.CALC) 

S2 

- 

(See  KBPR.CALC) 

X2 

= 

Axial  distance  at  point  2 

Z13 

(See  KBPR.CALC) 

KBPR.CPB 

C13 

= 

(See  KBPR.CALC) 

IGES 

= 

1 

MU1 

(See  KBPR.CALC) 

MU3G 

« 

(See  KBPR.CALC) 

Ml 

= 

(See  KBPR.CALC) 

M3G 

= 

(See  KBPR.CALC) 

R13 

= 

(See  KBPR.CALC) 

SMU1 

= 

Sin  (MU1) 

ST1 

= 

Sin  (11) 

SI 

= 

(See  KBPR.CALC) 

S2 

= 

(See  KBPR.CALC) 

TP3AV 

s = 

Average  value  of  the  streamline  angle  at  3 

Z13 

(See  KBPR.CALC) 

KBPR.SEPRTE 

MINF 

a 

Mach  number  (free  stream) 

PRT 

as 

Test  pressure  ratio 

PSOPI 

a 

^sep^ox 

PSOPIA 

s 

Psep/P® 

PSURF 

s 

(See  KBPR.CNTROL) 

PTOPI 

- 

Po/P« 

XSURF 

* 

(See  KBPR.CNTROL) 

KBPR. SETUP 

DEL1NC 

New  starting  position  for  the  characteristic  which  ends 

at  the  end  of  a  section 

D1 

K 

Starting  position  of  characteristic  at  initial  profile 

EKJ 

X 

Used  to  increment  J 
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JLIHT 


JST 

M 

NU 

NU2 

R 

RTO 

R2 

T 

T2 

X 

XI 
X2 
X3 

ACC 

BETA 

BETA1 

BETA2 

DIFF 

DNU 

DNUB 

G1 

G2 

L 

Ml 

M2 

NU 

NUMAX 

NU1 

NU2 

THETA 


(See  KBPR.CNTROL) 

Subscript 

Mach  number 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

Used  to  determine  point  closest  to  end  of  the  region 
(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

X(I-l.l) 

(See  KBPR.CNTROL) 

X( 1+1,1) 

KBPR.PMTURN 

Accuracy  desired 
Defined  by  equation 

/k*  -  1 
SG2 

BETA1  -  BETA2  (used  to  determine  convergence) 
Function  statement 
DNU  evaluated  at  BETA2 
/(G-l) / (G+l) 

1/G1 

Iteration  number 
Mach  number  upstream 
Mach  number  downstream 

Prandtl-Meyer  turn  angle  (function  statement) 

Maximum  value  of  Prandtl-Meyer  turn  angle 
NU  evaluated  at  BETA1 
NU1  +  THETAR 

Isentropic  turning  angle  in  degrees 
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THETAR 


Isentropic  turning  angle  in  radians 


A 

BETA 

G 

PMANGL 

AB 

AIAC 

CEF 

DX 

I  AC 

K 

L 

R1 

SOLI 

S0L2 

T 

X 

XHIGH 
XL 
XL  2 
XU 
XU  2 

A 

B 

C 

D 

DD 


KBPR. PMANGL 

/  (G-l) / (G+l) 

/M2  -  1 

Mach  number  upstream 
Prandtl-Mever  turn  angle 

KBPR.S0TE2B 

Used  as  a  final  check  for  accuracy 

Defined  by  equation 

EI/EJ  (checks  number  of  solutions) 

Increment  of  X 

AIAC 

Exponent 

Determines  which  root  is  desired  (+1  =  upper  root 
-1  *  lower  root) 

XU  -  XL 

Value  of  function  at  XU 
Value  of  function  at  X 
SOLI  *S0L2 

Independent  variable 
Defined  by  equation 
Lower  limit  of  search 
X 

Upper  limit  of  search 
XHIGH 

KB PR. SHOCK 
Defined  by  equation 
-  DM*D 

Defined  by  equation 
1/TD 

Detachment  turn  angle 
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G  =  Ratio  of  specific  heats  (=  Cp/Cv) 

G1  =  (G-D/2 

G2  =  (G+l)/2 

MIS  =  Ml2 

P  =  Defined  by  equation 

PHI  =  1/3  cos"1  (PHIARG) 

PHIARG  =  Defined  by  equation 

P3  =  Constant 

TD  =  Tan  (DE) 

X  =  Defined  by  equation 

Y  =  Defined  by  equation 

KBPR.FSTEP4 

CI2  =  Integral  12  evaluated  from  -3  to  +3 

CJ2  =  Integral  J2  evaluated  from  -3  to  +3 

OR  =  Ration  of  Crocco  numbers  C3a/C2a 

CS  =  C3a2 

EGP  =  Defined  by  equation  (estimated  geometric  parameter) 

ETAM3D  =  Three-dimensional  value  of  Dm 

ETAJ  =  Value  of  n  at  separating  streamline 

FSTEP4  =  EGP  -  GP 

GP  =  Geometric  parameter 

VII  =  Integral  II  evaluated  from  nj  to  +3 

VJ1  =  Integral  J1  evaluated  from  Dj  to  +3 

KBPR.BASE5 

ACC5C6  =  Accuracy  requirement  for  base  pressure  solution 

ACST4  =  Accuracy  requirement  for  KBPR.S0TE2B  subroutine 

C3AS  =  C33  (C3.J  “  Crocco  number  at  location  3a  in  base 

pressure  model) 

DXLIN  =  Increment  of  X  used  in  KBPR. LINEAR 
G  «  Ratio  of  specific  heats 


IBOUND  *  (See  KBPR.CNTROL) 

IOPT  »  Printing  option  number 

MIA  =  (See  KBPR.CNTROL) 

M2A  *  Mach  number  at  position  2a  in  base  pressure  model 
M3AM2A  =  M3A/M2A 

N  =  Number  of  iterations  in  KBPR. LINEAR 

NASHF  *  Nash's  recompression  coefficient 

PBPI  =  Base  pressure  ratio,  P^/Pl 

PEST  =  Estimated  base  pressure  ratio 

PHIJ3  =  Velocity  ratio  along  j-streamline 
THETD  =  THET12  in  degrees 
THET12  =  (See  KBPR. FLARE) 

THET3A  -  (See  KBPR. FLARE) 

THET3D  =  THET3A  in  degrees 
TRB  =  (See  KBPR. All) 

KBPR.FBASE5 

ACST4  *  (See  KBPR.FBASE5) 

BS  »  (See  KBPR.CNTROL) 

CI2  «  Integral  12  evaluated  from  -3  to  +3 

CJ2  =  Integral  J2  evaluated  from  -3  to  +3 

CR  *  Crocco  number  ratio,  C3A/C2A 

CS  =  C3A2 

Cl  =  Crocco  number  at  1 

C2A  =»  Crocco  number  at  location  2a  in  base  pressure  model 
C2AS  =  C2A2 

C3A  =  Crocco  number  at  location  3a 

C3AC2A  -  C3A/C2A 
C3AS  -  C3A2 

C5  “  Defined  by  equation 

C6  »  Defined  by  equation 

DINC  -  (See  KBPR.CNTROL) 

EGP  -  (See  KBPR.FSTEP4) 


ETAJ3 

ETAM2D 

ETAM3D 

FP 

G 

GP 

Cl 

G2 

G3 

G4 

I BOUND 

10UT 

MESHPM 

MIA 

MIAS 

M2A 

M2AS 

M3AM2A 

M3A 

M3AS 

PHI 

PHIJ3 

PNASH 

P1P3 

P2P3 

P4P3 

ROBDY 

RBSRSF 

RSRB 

THET12 

THF.T3A 

TR 


ETAJ  at  location  3  (see  KBPR.FSTEP4) 

Value  of  nm  for  two-dimensional  flow 
(See  KBPR. FSTEP4) 

(See  KBPR.BASE5 

Ratio  of  specific  heats 

Geometric  parameter 

(G-l)/2 

(G+l)/2 

(G-D/G 

G/(G-1) 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

MIA2 

(See  KBPR. BASES) 

M2A2 

M3A/M2A 

Mach  number  at  location  3a 
M3A2 

Velocity  ratio 

Velocity  ratio  along  j-streamline 

Pressure  with  Nash's  recompression  coefficient 

P1/P3 

P2/P3 

P4/P3 

(See  KBPR.CNTROL) 

(See  KBPR.CNTROL) 

(See  KBPR. FLARE) 

(See  KBPR. FLARE) 

(See  KBPR. FLARE) 
l./TRB  (see  KBPR. All) 
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TRB  =  (Sec  KBPR.AI1) 

XBASE  =  (See  KBPR.CNTROL) 

KBPR.SEPRTE 

MINF  =  Mach  number  (free  stream) 

PRT  =  Test  pressure  ratio 

PSOPI  =  Psep/Pa> 

PSOPIA  =  Psep^“  at  ^ree  stream  Mach  number 

PSURF  =  (See  KBPR.CNTROL) 

PTOPI  -  P0/Poo 

XSURF  =*  (See  KBPR.CNTROL) 


APPENDIX  B 


OPERATING  INSTRUCTIONS 

To  operate  this  computer  program,  a  minimum  of  eleven  cards 
must  be  read  in.  The  maximum  number  of  cards  is  a  function  of 
the  type  of  Mach  number  profile  desired  as  input. 

1)  READ  BS,  XSTART ,  XFLARE ,  PHI,  RBSRSF,  RBODYA,  IOPTL 

BS  is  the  shape  of  the  afterbody  which  is  desired.  This 
program  has  a  capability  to  analyze  only  a  conical  flare  or 
boattail.  This  corresponds  to  BS  -  1.  Any  other  read-in 
will  result  in  a  termination  of  the  program.  It  should  be 
noted  that  the  term  BS  has  been  declared  to  be  an  integer 
number. 

XSTART  is  the  X-coordinate  which  locates  the  beginning  of  the 

characteristics  solution  over  the  body.  This  XSTART  position 
falls  on  the  cylindrical  part  of  the  body  of  revolution. 
XSTART  need  not  be  zero.  XSTART  is  read  in  the  machine  as  a 
real  number.  (See  Fig.  9) 

XFLARE  is  the  X-coordinate  which  locates  the  beginning  of  the 

afterbody.  (XFLARE  _>  XSTART).  This  again  is  a  real  number. 
For  XFLARE  *  XSTART,  no  characteristics  solution  is  done  over 
the  cylindrical  portion  of  the  body,  but  Instead,  the  charac¬ 
teristics  begin  immediately  before  the  afterbody  region. 

(See  Fig.  9) 

PHI  is  the  angle  of  the  flare  or  boattail  in  degrees.  This  is 


solution. 


(C)  1,1.00,  1.00,  5.00,  2.00,  1.00,  2 


This  corresponds  to  a  body  with  a  5°  flare  whose  afterbody 
length  is  2.0  radii. 

2)  READ  G,  M(l) ,  OMACH 

G  is  the  ratio  of  specific  heats  (■  Cp/C^).  This  is  a  real 
number.  For  air,  G  =  1.40. 

M (1)  Is  the  surface  Mach  number  at  the  point  at  which  the 
characteristics  solution  begins. 

OMACH  is  the  free  stream  Mach  number.  Actually  this  may  be 
any  reference  Mach  number  and,  for  example,  may  refer  to  the 
Mach  number  ahead  of  the  shock  for  a  projectile. 

SAMPLE  CARD: 

1.40,  2.50,  3.00 

3)  READ  IREADl ,  IREAD2 ,  10UT 

I READ 1  determines  how  the  initial  Mach  number  profile  is  to 
be  set  up  for  analysis.  IREADl,  an  integer  number,  corresponds 
to  the  position,  XSTART. 

IREAD2  determines  how  the  profile  will  be  transferred  from 
the  cylindrical  body  to  the  afterbody.  This  corresponds  to 
the  position  XFLARE. 

I OUT  determines  how  the  profile  will  be  transferred  from  the 
afterbody  region  to  the  base  region.  This  corresponds  to  the 
position  XBASE. 
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H 


General  Coimnents  on  I READ 1 ,  IREAD2,  TOUT 
These  are  all  Integer  numbers  which  are  defined  for  0,  1,  2. 
If  IREAD1  »  0,  then  a  uniform  profile  is  assumed  at  XSTART.  If 
IREAD1  «  l,  then  a  profile  along  a  line  perpendicular  to  the  body 
is  read  in.  If  IREAD  1  m  2,  then  a  profile  is  read  in  which  is 
along  a  characteristic. 

It  must  be  noted  that  if  XFLARE  ■  XSTART,  then  IREAD1  = 
IREAD2.  Similarly,  if  XSTART  -  XFLARE  =  XBASE,  then  I RE ADI  » 
IREAD2  -  I0UT. 

SAMPLE  CARDS: 

(A)  0,  1,  2 

Uniform  profile  at  XSTART;  transfer  from  the  cylindrical 
body  to  the  afterbody  along  a  line  which  is  nearly  perpend¬ 
icular  to  the  body;  transfer  to  base  region  along  a  charac¬ 
teristic  line. 

(B)  0,  0,  0, 

This  card  would  represent  a  base  pressure  so)ution  with  a 
uniform  profile  read-in. 

CAUTION ;  Only  an  irrotational  profile  may  be  read  into  the  pro¬ 
gram.  A  rotational  Mach  number  profile  (e.g.,  a  viscous 
profile)  cannot  be  used  correctly  since  the  characteristics 
solution  requires  an  irrotational  flow. 

4)  READ  MESHPM,  CDEL1B,  CDEL1P 
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MESHPM  Is  an  integer  number  which  determines  the  number  of 


individual  discrete  turns  into  which  the  boattail  angle  is 
divided.  Usually,  MESHPM  »  2  is  sufficient.  This  is  also 
used  in  the  base  pressure  solution  for  the  same  reason. 

CDEL1B  is  a  real  number  which  sets  the  Increment  of  the 
first  point  off  the  surface  along  a  characteristic.  CDEL1B 
is  used  in  the  base  region. 

CDEL1P  is  again  a  real  number  which  sets  the  increment  of 
the  first  point  off  the  surface  along  a  characteristic  on 
the  body. 

General  Comments: 

CDEL1B  *  0.01  and  CDEL1P  “  0.10  is  sufficient  for  good  accuracy. 

5)  READ  JLIMT 

JLIMT  is  the  number  of  cards  which  will  be  read  into  the 
program  which  describes  the  initial  Mach  number  profile.  This 
card  is  not  read  if  IREAD1  *  0. 

6)  READ  X2(I),  R2  (I),  MACH, T2  (1) 

X2  (I)  is  the  X-coordinate  of  various  points  in  the  initial 
profile . 

R2  (I)  is  the  R-coordinate  of  various  points  in  the  initial 
profile . 

MACH  is  the  Mach  number  at  the  point  X(I) ,  R(I). 

T2  (I)  is  the  streamline  angle  (in  degrees)  at  the  point 
X(I),  R(I) . 
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General  Comments 


There  must  be  a  total  of  JLIMT  cards  read  Into  the  program 
with  each  card  containing  the  above  quantities  for  a  single 
point.  The  first  card  must  describe  conditions  on  the  sur¬ 
face  and  subsequent  cards  proceed  outward  from  the  body. 

There  must  be  at  least  two  cards  read  into  the  program,  and 
the  maximum  number  of  cards  is  60. 

7)  READ  31 

8)  READ  32 

9)  READ  33 

These  three  read  statements  are  run  identification  statements. 
These  may  be  punched  with  alpha-numeric  characters  not  ex¬ 
ceeding  55  spaces.  These  three  data  cards  will  be  printed  at 
the  beginning  of  the  base  pressure  section. 

10)  READ  DXLIN,  ACST4,  ACC5C6,  NASHF,  I0PT 

DXLIN  is  a  real  number  which  is  used  in  the  LINEAR  subroutine. 
A  function  is  evaluated  at  some  value  X  and  once  again  at  the 
value  X-DXLIN.  A  typical  value  for  DXLIN  is  0.005.  However, 
larger  values  may  be  used  if  so  desired. 

ACST4  is  an  accuracy  requirement  which  is  used  in  the  S0TE2B 
subroutine.  This  determines  to  what  accuracy  a  function  is  to 
be  solved.  A  reasonable  value  of  ACST4  is  0.005. 
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ACCSC6  is  an  accuracy  requirement  used  in  determing  the  base 


pressure  ratio,  P./P, .  A  typical  value  for  ACCSC6  is  0.001. 
b  1 

NASHF  is  Nash's  recompression  coefficient  which  has  been 
assumed  to  be  1.0. 

IOPT  is  a  printing  option  number.  If  IOPT  «  2  then  the  base 
pressure  parameters  are  printed  out.  If  IOPT  +  2,  then  print¬ 
ing  is  suppressed. 

11)  READ  IBOUND 

I BOUND  is  an  integer  number  which  determines  which  type  of 
base  pressure  solution  is  desired.  IBOUND  *  1  for  a  conetail 
solution  or  IBOUND  *  2  for  a  constant  pressure  solution. 

12)  READ  PEST,  TRB 

PEST  is  the  estimated  base  pressure  ratio,  P^/P^.  For  most 
reasonable  guesses  at  the  base  pressure  at  moderate  flare  or 
boattail  angles,  the  program  does  converge.  However,  for 
large  angles,  coupled  with  a  poor  guess,  the  program  may 
terminate.  It  is  recommended  that  a  guess  be  made  which  is 
close  to  the  assumed  answer,  since  the  calculation  procedure 
is  relatively  lengthy. 

TRB  is  the  temperature  ratio,  T./T  ,  where  T,  is  the  base 
-  b  oa  b 

temperature  and  is  the  free  stream  stagnation  temperature. 
General  Notes; 

This  read  statement  is  also  used  as  a  control  statement.  If 
PEST  <  0,  the  program  starts  over  and  a  new  case  is  run. 
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However,  is  TRB  <  0,  then  a  new  base  pressure  solution  may  be 
obtained  in  which  the  temperature  ratio  may  be  varied.  This, 
of  course,  implies  that  the  flow  along  the  body  remains  the 


same . 

For  example  if  Just  one  case  is  to  be  run,  the  last  two  cards 
will  be: 

0.820,  1.00 

-1.0,  1.0 

However,  if  two  temperature  ratios  are  desired,  then  an  addi¬ 
tional  card  is  used: 

0.820,  1.00 
0.900,  0.80 
-1.0,  1.0 


I 
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APPENDIX  C 


BOATTAI L/ FLARE- BASE5 
COMPUTER  PROGRAM  INPUT  AND  OUTPUT 
FOR  THE  ABERDEEN  PROVING 
GROUND  UNIVAC  1108 
(SAMPLE  CASES  INCLUDED) . 


SAMPLE  INPUT-OUTPUT  FOR 

UNIFORM  APPROACH  FLOW  (M  =3.0) 

— - -  c 

AND  A  ONE  CALIBER  7°  CONICAL 

BOATTAIL  WITH  y  =  1.4  AND  Tu/T  =  1.0 

b  oa 


1VT8*K3PR  i.  11 .  BPDT 

1  1,0. no. o.oo,  j.oon,  i  .000.2 

2  1.40,5.00,5.00 


5  O.J.l 

4  :, o.oi.o.io 

5  30C  BOATTAIL  3ASE  PRESSURE 

S  M-3.00,  PHIa-T .00  DEG 

■>  6  CAL  MOOEL.  1  CAL  BOATTAIL 

3  0.005,0.005,0.001,1.0,2 


9 

1 

10 

0.65,1.00 

11 

-  L . 

Q  ,1.0 

IMAP. S 

MAP28R1 

PL71-3 

03/25/7?  14: 

l. 

IN 

KBPR. MAIN 

2 . 

IN 

KBPR.  FLARE 

5 . 

IN 

kbpr.lchar 

4. 

IN 

K8PR.SIMR 

S. 

IN 

KBPR.DINCR 

6. 

IN 

KBPR. AXIL 

IN 

KBPR. A 1 2L 

8. 

IN 

KBPR.. -UK. 

9. 

IN 

KBPR..U2L 

10. 

IN 

KBPR. 11 

11. 

IN 

KBPR. 12 

12 . 

IN 

KBPR.J1 

13. 

IN 

KBPR.J2 

14. 

IN 

KBPR. BODY 

15. 

IN 

KBPR . TAB 

16. 

IN 

KBPR. ERF 

17. 

IN 

KBPR. PRES 

13. 

r.v 

KBPR.RS8F 

19. 

IN 

KBPR. SIGMA 

20. 

IN 

KBPR. LINEAR 

21. 

IN 

KBPR. CALC 

IN 

KBPR. SURF 

25. 

IN 

KBPR.CPB 

24. 

IN 

K8PR.SEPRTE 

25. 

IS 

KBPR. SETUP 

26. 

IN 

KBPR . PMTURN 

2". 

IN 

KBPR.  PMANGL 

23. 

IN 

K8PR.SOTE28 

29- 

IN 

KBPR. SHOCK 

30. 

IN 

KBPR.FSTEP4 

31. 

IN 

KBPR. BASES 

32. 

IN 

KBPR.FBASE5 
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M(  INFINITY") 

PHI  *  -7.00 

*5.00 

0EGR£ES 

PriNFINITYV'PT  *  ■ 

UNIFORM  PROFILE  WITH  M  *  3.0 

X 

M 

P/PT 

CP 

.000 

3.000 

.0272 

.0000 

.001 

3. 188 

.0206 

-.0387 

.001 

3.391 

.0153 

- .0694 

.001 

3.591 

.0153 

- .  0694 

.005 

3.391 

.01S3 

-.0694 

.109 

3.387 

.0154 

-.0689 

.  162 

3.383 

.0155 

-  .  0684 

.216 

3.379 

.0156 

-.06-8 

.2'0 

3.3'5 

.0157 

-.0673 

.  523 

3.371 

.0158 

-  .  0668 

.377 

3.367 

.0159 

-  .0665 

.430 

3.363 

.0159 

-  .0658 

.484 

3.359 

.0160 

-.0652 

.538 

3.356 

.0161 

- .0647 

.591 

3 . 352 

.0162 

-.0642 

.  64S 

3.548 

.0163 

-.0637 

.699 

3.344 

.0164 

-.0631 

.752 

3.340  ' 

.0165 

- . 0626 

.806 

5 . 536 

.0166 

-.0621 

.859 

3.333 

.0167 

-  .0615 

.913 

3.329 

.0168 

-  .0610 

.967 

3.525 

.0169 

- . 0605 

1.020 

3.321 

.0169 

-.0599 

1.074 

3.317 

.0170 

-.0594 

1.193 

3.314 

.0171 

-.0588 

1.312 

3.305 

0173 

-.0576 

1.432 

3.297 

.0176 

-.0564 

1.S51 

3.289 

.0178 

-.0552 

1.670 

3.281 

.0180 

-.0539 

1.790 

3.272 

.0182 

-.0527 

1.909 

3.264 

.0184 

-.0514 

2.000 

3.256 

.0186 

-.0501 

2.000 

3.250 

.0189 

-.0491 
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SAMPLE  INPUT-OUTPUT  FOR  THE 


NON-UNIFORM  APPROACH  FLOW 

OBTAINED  FROM  BRL  METHOD  OF 

CHARACTERISTICS  PROGRAM  AND 

A  ONE  CALIBER  7°  CONICAL  BOATTAIL 

WITH  y  =  1.4  AND  T./T  =1.0. 

b  oa 


WTR*KflPR( 1 ) . 3PDAT 

1  1 ,0 .00 .0 .00 . -7 .00,2.000,1 .00.2 

2  1.40.3.135,3.00 

3  1,1.2 

4  2,0.01,0.10 

5  22 

6  0.000.1.000,3.050,0.000 

7  0. 000, 1.000, 3. OSO, 0.000 

8  0.000, 1.200, 3. OSO, 0.000 

9  0.000.1.300,3.070,0.000 

10  0.000,1.400,3.100,0.000 

11  0.000,1.500,3.080,0.000 

12  0.000,1.500,3.050,0.000 

13  0.000,1.700,5.000,0.000 

14  0.000.1.800,2.960,0.000 

15  0.000,1.900,2.920,0.000 

16  0.000,2.000,2.890,0.000 

17  0.000.2.100,2.860,0.000 

18  0.000,2. 200, 2. 8S0, 0.000 

19  0.000,2.300,2.830,0.000 

20  0.000,2.400,2.820,0.000 

21  0.000,2.500,2.860,0.000 

22  0.000,2.600,2.970,0.000 

23  0.000,2.700,3.000,0.000 

24  0.000,2.800,3.000,0.000 

25  0.000,2.900,3.000,0.000 

26  0.000,3.000,3.000,0.000 

27  0.000,3. 100,3.000,0.000 

28  SOC-BOATTAIL  BASE  PRESSURE 

29  M-3.00.  PHI--7.00  DEG 

30  6  CAL  MODEL  -  l  CAL  BOATTAIL 

31  0.005,0.005,0.001,1.0,2 


52  1 

53  0.750,1.00 

34  -1.0, 1.0 
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APPENDIX  D 


NOTATION: 


COMPUTER  PROGRAM 
FLOW  CHARTS 


1.  Main  Program  or  subprogram  |  j 

2.  "COMMON"  block  Q  !  information  is 
transferred  between  subf5rograms  FI  and  F2 
through  COMMON  block  Cl. 


3.  Subprogram  CALL'S 


B  CALL’S  C  as  explicitly  Instructed  in  A. 
That  is,  C  is  an  EXTERNAL  function  to  A's 
CALL  to  fi. 


The  order  of  execution  generally  proceeds  from  top  to  bottom  and/or  from 
left  to  right.  Internal  functions  such  as  SQRT ,  T  AN ,  etc.  are  not  shown. 
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DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


No.  of 
Copies 


Organization 


12  Commander 
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